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A  Model  for  High  Frequency  Radar  Auroral  Clutter 


I.  THE  RADIO  AURORA 

1.1  Introduction 

The  so  called  "radio-aurora"  or  "radar-aurora"  is  a  phenomenon  closely  re¬ 
lated  to  the  visible  aurora,  but  differing  from  the  latter  in  several  important  re¬ 
spects.  The  visible  aurora  is  caused  by  optical  emissions  stimulated  by  collisional 
excitation  of  atmospheric  atoms  by  energetic  charged  particles  of  magnetospherie 
origin.  The  radio  aurora,  on  the  other  hand,  is  the  result  of  radio  wave  scattering 
from  inhomogeneities  in  the  ionospheric  plasma.  It  happens  that  the  general  region 
of  the  visible  aurora  is  also  one  of  intense  small  scale  ionospheric  irregularities, 
so  that  there  are  similarities  in  the  spatial  and  temporal  characteristics  exhibited 
bv  the  visible  and  radio  auroras.  However,  since  these  phenomena  are  generated 
bv  different  physical  processes,  it  may  reasonably  be  expected  that  they  would 
exhibit  some  different  properties. 

The  morphology  of  the  visible  aurora  has  been  extensively  studied,  since  the 
necessary  optical  techniques  have  been  available  for  over  a  century.  It  has  become 
well  accepted  that  the  visible  aurora  is  confined  to  a  narrow  region  (a  "ring"  or 
"belt")  surrounding  the  geomagnetic  poles.  With  the  relatively  recent  advent  of 
sensitive  photometers  and  satellite  optical  scanners,  a  hitherto  unsuspected  optical 
component  of  the  aurora  was  discovered  which  is  below  the  intensity  level  for  normal 
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optical  or  photographic  viewing.  This  diffuse,  sub-visual  aurora  occupies  a  greater 
region  of  space  than  the  bright  visible  aurora,  and  appears  to  be  more  stable  in 
time.  Since  the  charged  particle  precipitation  responsible  for  the  visible  and  sub¬ 
visual  aurora  is  strongly  affected  by  the  geomagnetic  field,  these  phenomena  are  * 

best  ordered  in  a  system  of  coordinates  based  on  the  structure  of  the  earth's  mag¬ 
netic  field.  Various  geomagnetic  coordinate  systems  have  been  proposed  and  used 
to  order  charged  particle  precipitation  observations.  The  one  currently  most  used  ’ 

is  that  of  Hakura,  *  known  as  the  "corrected  geomagnetic"  coordinate  system  (CGS), 
which  is  discussed  in  Section  1. 4. 

The  radio  aurora  is  also  better  ordered  in  the  CGS  than  in  the  conventional 
geographic  coordinate  system,  but  it  does  not  exhibit  precisely  the  same  properties, 
in  this  system,  as  does  the  visible  aurora.  The  ionospheric  inhomogeneities  re¬ 
sponsible  for  radio  wave  scattering  are  caused  by  plasma  instabilities  which  are 
driven  by  ionospheric  currents.  These  currents  in  turn  flow  in  response  to  im¬ 
pressed  electric  fields  of  magnetospheric  origin,  and  their  intensity  and  space- 
time  distribution  reflect  both  the  properties  of  these  impressed  electric  fields  and 
of  the  conductivity  of  the  ionosphere.  The  conductivity  depends,  among  other  things, 
on  the  degree  of  ionization  of  the  atmosphere,  and  to  the  extent  that  this  ionization 
is  caused  by  particle  precipitation,  this,  too,  is  best  ordered  in  the  CGS.  However, 
solar  ultra-violet  radiation  also  ionizes  the  atmosphere  to  an  extent  dependent  on 
the  elevation  angle  of  the  sun,  so  there  is  a  component  of  ionospheric  conductivity 
which  is  ordered  in  the  standard  geographic  coordinate  system.  This  complexity 
is  partly  responsible  for  the  different  properties  of  visible  and  radio  auroras. 

Intense  ionospheric  currents  do  not  necessarily  flow  in  the  same  region  of  the 
ionosphere  as  that  in  which  intense  particle  precipitation  occurs.  Thus,  although 
information  relating  to  charged  particle  precipitation  and  optical  emissions  is 
available  in  profusion,  care  needs  to  be  exercised  in  using  this  data  to  deduce 
properties  of  the  radio  aurora. 

In  this  section,  the  relationship  between  small  scale  ionospheric  irregularities 
and  radio  auroral  clutter  is  briefly  reviewed.  A  review  of  magnetosphere-ionos¬ 
phere  interaction  follows,  aimed  at  placing  the  subject  of  the  high  latitude  ionos¬ 
phere  in  proper  geophysical  perspective.  The  various  coordinate  systems  employed 
in  modelling  the  high  latitude  ionosphere,  auroral  phenomena,  and  radar  clutter 
are  discussed  and  the  concept  of  the  "auroral  oval"  is  reviewed.  Finally  there  is 
a  discussion  of  the  regular  or  predictable  features  of  the  high  latitude  ionosphere. 

The  model  of  auroral  radar  clutter,  which  will  subsequently  be  developed,  is  based 
on  the  above  concepts. 

1.  Hakura,  Y.  (1965)  Tables  and  maps  of  geomagnetic  coordinates  corrected  by 

the  higher  order  spherical  harmonic  terms,  Kept.  Ionosph.  Space  Res.  » 

(Japan),  19:121. 
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1.2  Ionospheric  Irregularities  and  Clutter 

Radar  clutter,  induced  by  the  aurora,  is  simply  the  result  of  radio  wave  back- 
scatter  from  ionospheric  irregularities  in  the  vicinity  of  the  aurora.  These  ir¬ 
regularities  are  strongly  anisotropic  in  form,  due  to  the  anisotropy  of  the  ionos¬ 
pheric  plasma  conductivity  brought  about  by  the  presence  of  the  geomagnetic  field. 
The  irregularities  may  be  visualized  as  thin  plasma  cylinders,  whose  axes  are 
aligned  with  the  local  geomagnetic  field  direction,  and  present  in  a  continuum  of 
scale  sizes.  Due  to  the  high  degree  of  anisotropy  in  the  shape  of  the  irregularities, 
their  radio  backscatter  cross  section  is  highly  directional,  that  is,  it  is  a  strong 
function  of  the  angle  of  incidence  of  the  radio  propagation  vector  on  the  geomagnetic 
field  direction. 

The  two  prerequisites  for  the  occurrence  of  auroral  clutter  are: 

(1)  A  region,  in  the  ionosphere,  of  small  scale  irregularities 
must  be  illuminated  by  the  radar, 

(2)  The  radio  propagation  vector  must  be  approximately 
orthogonal  to  the  geomagnetic  field  vector  in  the  region 
of  ionospheric  irregularities. 

Since  the  spatial  and  temporal  distributions  of  ionospheric  irregularities  are  mod¬ 
erately  well  known,  the  incidence  of  auroral  clutter  (and  techniques  for  minimizing 
it)  reduce  largely  to  a  study  of  the  factors  influencing  the  magnetic  aspect  angle. 
These  factors  are  primarily  the  radio  frequency  and  the  ambient  ionospheric  struc¬ 
ture,  which  together  determine  the  amount  of  refraction  which  the  radar  beam 
experiences. 

Were  it  not  for  the  large  refraction  which  the  ionosphere  imposes  on  HF  radio 
waves,  the  regions  of  near-orthogonal  magnetic  aspect  angle  would  be  very  small, 
as  they  are  in  the  case  of  UHF  line-of-sight  propagation.  Even  if  the  HF  radar 
possessed  a  narrow  ("pencil")  antenna  beam,  the  assessment  of  effects  of  auroral 
clutter  and  the  means  of  its  mitigation  would  not  present  major  problems.  How¬ 
ever,  the  fact  that  the  radar  beam  is  usually  very  wide  in  the  elevation  plane  adds 
considerable  complexity  to  the  problem  since  the  ionospheric  refraction  then  per¬ 
mits  the  orthogonality  criterion  to  be  met  over  a  very  large  range  extent.  In  fact, 
when  the  possibility  of  multiple  ionospheric  hops  is  considered,  the  range  extent 
over  which  the  orthogonality  criterion  can  be  satisfied  becomes  virtually  unlimited. 
This  does  not  imply  that  the  amplitude  of  the  observed  auroral  clutter  will  be  uni¬ 
form  in  range;  several  factors,  principally  radio  wave  absorption,  contribute  to 
pronounced  variations  of  clutter  amplitude  with  range.  In  particular,  in  the  ease 
of  multiple  hop  propagation,  there  may  be  a  large  amplitude  difference  between 
clutter  signals  observed  via  one  and  two  or  more  hops,  especially  in  daytime  when 
the  absorption  in  the  lower  ionosphere  is  strongest. 


1.3  Magiiftosplierit’  lonospliere  Interaction 

An  understanding  of  the  special  features  of  the  high-latitude  ionosphere  requires 
an  appreciation  of  the  relationship  between  ionospheric  and  magnetospheric  proper¬ 
ties.  The  earth's  magnetic  environment  (loosely  termed  the  "magnetosphere")  is 
in  many  ways,  analogous  to  a  plasma  bottle  in  laboratory  magneto-plasma  physics. 
Figure  1  shows  a  schematic  representation  of  the  magnetosphere.  The  solar  wind- 
a  stream  of  charged  particles,  primarily  electrons  and  protons,  emanating  from 
the  sun— impinges  on  the  earth's  magnetic  field.  A  complex  interaction  takes  place 
in  which  a  standing  shock  wave  is  formed  at  a  geocentric  distance  of  about  15  earth 
radii  (R  ),  and  the  magnetic  field  lines  behind  the  shock  are  swept  back  to  form  a 
"tail"  in  anti-solar  direction.  The  region  immediately  behind  the  shock  contains 
hot  plasma,  moving  at  a  velocity  less  than  the  local  speed  of  sound;  this  plasma  has 
a  narrow  region  (cleft)  of  direct  access  to  the  dayside  of  the  ionosphere,  between 
the  first  swept-back  field  line  and  the  last  "unbroken"  field  line.  At  the  center  of 
the  tail  region  there  exists  a  narrow  sheet  of  plasma  which  is  in  constant  convective 
motion  and  in  which  some  spatial  energy  differentiation  is  maintained.  This 
"plasma-sheet"  is  also  directly  connected  to  the  high-latitude  ionosphere,  on  the 
nightside  of  the  earth  and  its  plasma  is  subject  to  acceleration  by  sporadic  elec¬ 
tric  fields  which  are  induced  across  the  tail,  in  the  Fast-West  direction,  by  the 
solar  wind.  Figure  2  shows  an  equatorial  cross  section  of  the  magnetosphere  in 
which  the  plasma  convection  pattern  is  illustrated.  In  both  Figures  1  and  2,  a 
region  known  as  the  "plasmasphe re"  is  depicted,  extending  to  a  geocentric  distance 
of  about  4  R  .  The  plasmasphere  boundary  is  the  surface  separating  plasma  which 
is  co-rotating  with  the  earth  from  that  which  is  being  converted  under  the  influence 
of  the  solar  wind.  For  many  purposes  the  demarcation  between  the  high-latitude 
ionosphere  and  mid-latitude  ionosphere  may  be  defined  as  the  intersection  in  the 
ionosphere  of  the  magnetic  field  lines  which  bound  the  plasmasphere. 

1.4  Coordinate  System* 

A  complicating  factor  in  studying  the  high-latitude  ionosphere  and,  specifically 
in  modelling,  it  is  the  fact  that  major  parameters  are  ordered  in  two  different  coor¬ 
dinate  systems.  In  computing  the  combined  effects  of  these  parameters  it  is  always 
necessary  to  keep  track  of  their  position  in  each  system.  In  particular,  when  com¬ 
puting  radar- related  quantities  such  as  clutter  intensity  and  Doppler  ship‘  of  the 
radar  signal,  it  is  of  course  necessary  to  consider  a  system  of  coordinates,  defined 
by  range  and  azimuth,  centered  on  the  radar.  The  temperate  zone  ionosphere  is 
largely  ordered  in  the  conventional  geodetic  system  of  coordinates,  defined  by  geo¬ 
graphic  latitude  and  longitude,  and  local  time.  This  follows  .rom  the  fact  that  solar 
illumination  is  the  dominant  contributing  factor  to  its  production  and  consequently 


10 


While  it  was  well-known,  at  an  early  date,  that  high-latitude  ionospheric 
properties,  especially  the  aurora,  were  best  organized  and  represented  in  geomag¬ 
netic  coordinates,  substantial  discrepancies  and  anomalies  were  observed  when  the 
earth's  magnetic  field  was  approximated  by  a  dipole.  A  gradual  understanding  of 
the  importance  of  higher  order  terms  in  the  spherical  harmonic  representation  of 
the  true  earth's  magnetic  field  culminated  in  the  formulation  of  the  "corrected 
geomagnetic  (CG)  coordinate  system"  — abbreviated  in  this  work  as  C’GS.  The 
charged  particles  which  generate  auroral  and  ionospheric  effects  are  strongly  in¬ 
fluenced  by  magnetic  field  structure  at  large  distances  from  the  earth;  this  results 
in  a  magnetic  latitudinal  dependence  of  particle  precipitation,  together  with  a 
temporal  variation  introduced  by  the  earth's  rotation.  The  coordinate  system  best 
adapted  to  ordering  charged  particle  precipitation  data  (as  well  as  other  related 
effects)  is  defined  by  corrected  geomagnetic  latitudes  U'Gl.)  and  corrected  geomag¬ 
netic  time  (CUT).  For  some  approximate  purpose's,  a  reasonable  degree  of  accu¬ 
racy  is  obtained  by  substituting  local  (solar)  time  for  corrected  geomagnetic  time. 

The  transformation  between  the  geographic  coordinate  system  and  the  ('(IS  is 

•> 

accomplished  by  a  numerical  mapping  procedure,  based  on  the  work  of  (iustaffson.  “ 

I.S  The  Auroral  Oval 

A  major  advance  was  achieved  in  understanding  the  morphology  oftho  high-lati¬ 
tude  ionosphere  with  the  discovery  of  the  auroral  oval.  Prior  to  this  discovery,  it 
had  been  known  that  auroral  activity  was  confined  to  a  fairly  narrow  range  of  lati¬ 
tudes,  centered  on  the  geomagnetic  pole,  but  no  clear  concept  existed  of  the  coupling 
of  spatial  and  temporal  aspects  of  auroral  occurrence.  The  aurora  is  now  known  to 
be  caused  by  the  injection  of  energetic  charged  particles  into  the  earth's  atmos¬ 
phere,  and  its  spatial  extent  is  determined  by  the  projection  on  the  ionosphere  of 
the  magnetospheric  regions  of  plasma  having  direct  access  to  the  ionosphere  (see 
Figures  1  and  2).  Since  the  earth's  rotation  introduces  a' time  variability  factor, 
it  is  necessary  to  consider  the  auroral  region  as  having  both  latitudinal  and  temporal 
dependences.  The  auroral  oval  is  depicted  in  Figure  3,  showing  its  inner  and  outer 
boundaries  as  heavy  lines,  in  COS  coordinates  at  12  FT,  for  a  moderate  level  of 
geomagnetic  activity  (Q  =  3).  It  is  important  to  note  that  the  auroral  oval  concept 
is  a  statistical  one,  deduced  originally  from  photographic  data  and  describing  the 
visible  aurora.  Subsequent  work  has  shown  that,  with  minor  modifications,  the 
auroral  oval  is  a  valid  representation  of  some  other  effects,  in  particular  iono¬ 
spheric  structure,  and  that  the  concept  has  some  validity  on  an  instantaneous  basis 
as  well  as  a  statistical  one.  The  auroral  oval  is  a  valuable  device  for  ordering 

2.  Gustafsson,  G.  (1970)  A  revised  corrected  geomagnetic  coordinate  system, 

Arkiv  f.  Geofvsik,  5(No.  401:595. 

" - ■ 
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M,.p  of  Northern  Polar  Kegions  Showing  Location  of  the  Auroral  Oval  at  12  UT  for 
eomagnetic  Activity.  Locations  of  vertical  incidence  ionosondes  during  the  IGY  an 


and  describing  high  latitude  ionospheric  phenomena,  but  certain  limitations  must 

be  borne  in  mind.  Aside  from  the  fact  that  the  auroral  oval  is  a  descriptor  of  a 

certain  type  of  optical  emissions  (bright,  spatially  discrete,  and  limited  in  spectral 
3 

width),  Eather  has  pointed  out  that  the  concept  ignores  variations  in  intensity  of 
particle  precipitation  with  time.  Such  variations  around  the  oval  are  important  in 
assessing  the  impact  of  auroral  clutter  on  a  backscatter  radar  system. 

1.6  Regular  Ionospheric  Features  at  High  Latitude* 

Regular  features  of  the  high-latitude  ionosphere  may  be  defined  as  those  which 
conform  to  the  morphological  model  illustrated  in  Figures  1,  2,  and  3.  The  effects 
of  the  magnetosphere -ionosphere  coupling  are  superimposed  on  the  well  known 
direct  solar  ionization  of  the  atmosphere  which  is  manifested  globally.  The  latter, 
of  course,  is  best  represented  in  the  standard  geographic  coordinate  system  (geo¬ 
graphic  latitude  and  solar  time);  this  leads  to  complications  since  the  superposition 
of  the  two  controlling  influences  must  be  made  with  due  regard  to  the  relevant 
coordinate  systems.  It  is  also  necessary  to  consider  the  special  nature  of  the 
direct  solar  ionizing  influence  at  high  latitudes,  where  the  sun  can  be  above  or 
below  the  horizon  for  long  periods,  in  summer  and  winter  respectively.  The  prin¬ 
cipal  effects  of  this  feature  are; 

(a)  The  diurnal  and  seasonal  variations  of  ionospheric  properties  are 
strongly  coupled, 

(b)  In  winter,  in  the  absence  of  direct  solar  ultra  violet  (UV)  radiation, 
the  particle- related  effects  tend  to  predominate, 

(c)  In  summer,  the  polar  ionosphere  is  dominated  bv  solar  UV  ionizatio... 

While  the  roles  of  direct  solar  illumination  and  charged  particle  precipitation 

are  at  least  intuitively  apparent,  other  physical  factors  are  also  at  work  in  influenc¬ 
ing  the  structure  of  the  quiet  polar  ionosphere.  Perhaps  the  most  significant  of 
these  factors  are  the  chemical  composition  of  the  upper  atmosphere  and  the  coupling 
between  the  neutral  and  ionized  atmospheric  constituents.  The  chemical  composi¬ 
tion  of  the  high  atmosphere  in  the  polar  region  is  dominated  by  an  upwards  diffusive 
flow,  known  as  the  "polar  wind,"  which  results  in  a  great  reduction  of  the  lighter 
atmospheric  constituents  (hydrogen  and  helium)  relative  to  the  corresponding  com¬ 
position  at  mid-latitudes.  The  net  result  is  a  general  reduction  of  electron  density 
over  the  polar  cap,  together  with  an  increase  in  the  effective  plasma  scale  height, 
especially  at  high  altitudes.  The  viscous-like  coupling  between  movements  of 
neutral  and  ionized  atmospheric  constituents  gives  rise  to  ionospheric  effects  in¬ 
duced  by  neutral  winds  which  can  be  a  significant  factor  in  the  long  winter  nighttime 
when  ionization  production  is  low. 

T!  (lathe r,  R.  II.  (li'73)  The  auroral  oval  -  A  reevaluation.  Rev.  (leophvs.  Space 

Plus.  11:135. 

- * — 

Plasma  scale  height  is  a  parameter,  having  dimensions  of  length,  which  charac¬ 
terizes  the  exponential  decav  with  height  of  plasma  density. 


1.  6.  1  THF  IONOSPHERIC  TROUGH 

A  notable  feature  of  the  high-latitude  ionosphere  is  the  ionospheric  trough, 
which  is  a  region  of  strongly  depressed  electron  density,  occurring  at  F -region 
altitudes,  primarily  on  the  night-side  of  the  earth.  The  trough  is  aligned  approxi¬ 
mately  with  a  circle  at  a  constant  geomagnetic  latitude  of  about  60°,  The  equatorial 
edge  of  the  trough  coincides  approximately  with  the  locus  of  the  plasmapause,  and 
the  electron  density  there  increases  gradually  from  the  deepest  point  (~  60°  CGL) 
to  the  mid-latitude  value.  The  polar  edge  of  the  trough  is  much  steeper  than  the 
equatorial  edge  and  coincides  with  the  equatorial  edge  of  the  auroral  oval.  The  local 
time  dependence  of  the  electron  density  in  the  trough  region  shows  a  minimum  at 
about  0300,  and  onset  and  disappearance  at  variable  time,  corresponding  respective¬ 
ly  to  sunrise  and  sunset  at  the  trough  latitude.  Thus  the  trough  can  be  represented 
schematically  in  the  CG  coordinate  system  as  shown  in  Figure  4. 
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Figure  4.  Schematic 
Representation  of  the 
F -layer  Ionospheric  Trough 
and  Plasma  Ring,  in 
Corrected  Geomagnetic 
Coordinates,  Relative 
to  the  Auroral  Oval 


The  ionospheric  trough  is  a  highly  variable  feature  which  appears  to  be  very 
prevalent  in  the  topside  ionosphere  but  is  also  frequently  present  at  the  F-peak  and 

4 

below.  Spiro  et  al  have  presented  a  convincing  argument  explaining  the  trough  as 
the  result  of  the  complex  magnetospheric  convective  flow  near  the  plasmapause. 
Ionospheric  plasma  in  this  region  is  influenced  bv  the  earth's  rotation  and  bv  the 
cross -tail  electric  field  in  such  a  way  that  a  region  of  How  stagnation  exists  on  the 

4.  Spiro,  R.W.,  Heelis,  R.  A.  ,  and  Hanson,  W.  13.  ( 1 07 8 )  Ion  convection  and  the 
formation  of  the  mid-latitude  ionization  trough,  .1.  Geophys.  Res.  83;4255. 
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night  side  of  the  earth.  In  the  absence  of  ionization  sources,  normal  K- region 
recombination  acting  over  a  long  period  results  in  strong  reduction  in  the  degree 
of  ionization,  giving  rise  to  a  trough-shaped  plasma  depression.  The  trough  is 
observed  most  frequently  in  winter  and  in  the  early  morning,  and  its  intensity  and 
location  are  fairly  closely  correlated  with  the  level  of  geomagnetic  activity.  Thus, 
the  latitude  (CGL)  of  the  trough  electron  density  minimum  moves  equatorwards  by 
1.8  -  2.  0  degrees  per  unit  increase  in  Kp  (the  planetary  magnetic  index). 

1.6.2  AURORAL  E-LAYER 

In  the  vicinity  of  the  auroral  oval  there  appears  a  relatively  regular,  thick 
layer  of  ionization  at  E-region  altitudes,  produced  by  a  steady  flux  of  charged 
particles  in  the  kilo-electron  volt  (keV)  energy  range.  The  maximum  plasma  fre¬ 
quency  in  this  "auroral-E'Taver  is  typically  1.  5  -  2.  5  MHz.  On  the  day-side  of  the 
earth,  the  auroral-E  contribution  augments  the  solar  UV-produced  E-layer. 

1.6.3  AURORAL  SPORADIC-LAYER 

Highly  variable  fluxes  of  charged  particles  (electrons  and  protons)  produce 
auroral  sproadic-E,  which  is  mainly  confined  to  the  auroral  oval  region.  Plasma 
instabilities  induced  by  the  intense  "auroral  electrojet"  current  system  give  rise  to 
additional  small  scale  structure  in  the  E- region  within  and  immediately  adjacent  to 
the  auroral  oval.  Auroral  sporadic-E  is  capable  of  scattering  incident  radio  waves 
at  very  high  frequencies  (well  into  the  VHP  range)  and  create  serious  multipath 
problems  on  radio  links  in  the  auroral  region. 

1.6.4  F- REGION  IRREGULARITIES 

The  F-region  in  the  auroral  nighttime  ionosphere  very  often  displays  sub¬ 
stantially  lower  electron  densities  than  the  E-region,  and  is  virtually  always  irregu¬ 
lar  in  nature.  On  the  day-side  of  the  auroral  oval,  a  special  situation  is  observed 
virtually  permanently;  a  region  of  several  degrees  latitudinal  extent  containing 
irregular  F-region  ionization  of  relatively  high  intensity  is  observed,  which  is  the 
result  of  direct  access  of  solar  wind  plasma  into  the  ionosphere  through  the  mag- 
netospheric  cleft.  The  zone  surrounding  the  magnetic  pole  in  which  low  energy 
particle  precipitation  produces  excess  irregular  ionization  is  sometimes  referred 
to  as  the  "plasma  ring",  depicted  schematically  in  Figure  4.  The  region  interior 
to  the  plasma  ring,  in  the  immediate  vicinity  of  the  magnetic  pole,  is  one  of  depleted 
electron  density,  especially  during  the  long  winter  night;  this  region  is  sometimes 
referred  to  as  the  "polar  cavity.  " 
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R  MHO  PROPAGATION  CONSIDERATIONS 


An  essential  element  in  computing  the  incidence  and  intensity  of  auroral  clutter 
on  a  111-'  radar  is  the  computation  of  the  radio  ray  trajectories  in  the  ionosphere. 

This  involves  the  construction  of  a  suitable  numerical  model  for  the  ionosphere, 
appropriate  to  the  geophysical  conditions  of  interest,  and  subsequent  numerical  ray 
tracing  through  this  model.  In  order  to  compute  the  intensity  of  auroral  backscatter. 
a  model  for  the  ionospheric  irregularity  cross  section  is  necessary,  and  an  esti¬ 
mate  of  the  absorption  of  the  radio  waves  in  the  ionosphere  must  be  made.  Because 
of  the  aspect  sensitivity  of  auroral  backscatter,  a  suitable  global  model  for  the  geo¬ 
magnetic  aspect  angle  must  be  computed  at  all  points  along  the  ray  trajectory. 

2.1  lonoiiptierk’  Model* 

The  model  of  the  ionosphere,  chosen  in  this  study  to  represent  the  high-latitude 
ionosphere,  is  one  developed  at  Home  Air  Development  Center  and  its  predecessor 
organization— AFC  HI..  This  model  is  described  by  Hush  et  al5  and  Miller  and 

g 

llibbs.  The  dominant  component  of  any  ionospheric  model  intended  for  ray  tracing, 

is  the  f  K,,  model.  In  the  HADC  model,  this  component  is  essentially  that  of  the 
°  -  7 

1TS-78  model  with  empirical  modifications  to  incorporate  magnetic  activity  de¬ 
pendence  and  the  ionospheric  trough  in  the  CCS.  The  model  provides  a  three- 
dimensional  median  numerical  ionospheric  specification  as  a  function  of: 

(1)  Sunspot  number  (SSN). 

(2)  Planetary  magnetic  index  (Kp>, 

Cl)  Universal  lime  <UT). 

(4)  Season  (that  is.  month). 

(5)  l  atitude  and  longitude. 

2.2  Numerical  R»y -Tricing 

Huy  tracing  is  the  technique  of  computing  the  trajectories  of  radio  waves  through 

g 

a  modelled  ionosphere  by  solution  of  the  so-called  llaselgrove  equations.  This  is 
accomplished  numerically  by  use  of  a  computer  program  originally  developed  by 
.tones  and  later  modified  extensively  at  Home  Air  Development  Center.  ** 

Among  the  modifications  made  are  the  inclusion  of  sporadic-E,  the  input  of  arbi¬ 
trary  ionospheric  specifications,  the  calculation  of  deviative  absorption  and  a  vari¬ 
ety  of  output  options.  For  the  purposes  of  the  present  study,  the  modification  of 
prime  interest  is  the  incorporation  of  an  accurate  geomagnetic  field  model,  and  the 
continuous  computation,  along  the  ray  trajectory,  of  the  angle  between  the  propaga¬ 
tion  vector  and  the  geomagnetic  field  vector  (magnetic  aspect  angle).  The  geomag¬ 
netic  field  model  will  be  described  in  Section  2,  5. 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  deferences,  page  87.) 
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The  ray  tracing  program  permits  computation  of  the  ray  trajectory  over  an 
arbitrary  number  of  hops  for  a  given  azimuth  and  elevation  angle  from  any  chosen 
location.  In  the  present  study,  the  maximum  number  of  hops  was  taken  as  two,  and 
rays  were  computed  in  an  azimuthal  range  of  -60°  to  +60°,  and  for  elevations  of 
1°  to  21°. 

2.3  Absorption 

The  calculation  of  ionospheric  absorption  in  the  polar  regions  is  complicated 
by  the  presence  of  two  different  absorption  components: 

(1)  "Normal"  D- region  absorption,  which  is  primarily  a  function 
of  the  intensity  of  the  solar  ionizing  radiation  (I.vman-alpha) 
and  is  therefore,  dependent  on  the  solar  zenith  angle, 

(2)  "Auroral"  absorption,  caused  by  ionization  of  the  lower 
ionosphere  by  energetic  charged  particle  precipitation. 

The  "normal"  component  of  absorption  is  relatively  deterministic  and  exhibits 
a  large  seasonal  component  which  "seasonal  anomaly"  is  apparently  due  to  seasonal 

changes  in  the  molecular  composition  of  the  lower  ionosphere,  and  is  also  latitude- 

1 2 

dependent.  A  model  for  the  normal  D-  region  absorption  has  been  accepted  as  an 
international  standard  by  the  CCIR,  and  is  used  for  computation  in  this  study.  A 
computer  program,  based  on  the  Bradley ~G eo rge ^  algorithm,  was  incorporated 
into  the  ray  trajectory  calculations. 

The  auroral  component  of  absorption  is  more  complicated,  being  essentially  a 

random  process  but  possessing  certain  deterministic  features.  Among  the  many 

workers  who  have  studied  auroral  absorption,  Hartz  et  al13  and  Hargreaves  and 
14 

Cowley  have  published  quantitative  data  suitable  for  computational  purposes. 

15 

More  recently  Foppiano  has  published  a  numerical  method  of  predicted  auroral 

1 6 

absorption  of  HF  skywaves.  Elkins  has  summarized  the  quantitative  features  of 
auroral  absorption  events  in  a  manner  suitable  for  computation  for  a  typical  case. 
Auroral  absorption  occurs  during  so-called  "auroral  sub-storms,"  and  is  confined 
to  a  belt  of  latitude  approximately  5°  -  10°  in  extent,  located  slightly  equatorwards 

12.  Bradley,  P.  L. ,  and  George,  P.  A.  (1974)  A  new  method  of  predicting  the 

ionospheric  absorption  of  high  frequency  waves  at  oblique  incidence. 
Telecomm.  J.  41:307. 

-  .w. 

13.  Hartz,  T.  R. ,  Montbriand,  L.  E.  ,  and  Vogan,  E.  1..  (1963)  A  study  of  auroral 

absorption  at  30  Mc/s,  Can  J.  Phys,  41:581. 

14.  Hargreaves,  J.K. ,  and  Cowley,  F.C.  (1967)  Studies  of  auroral  radio  absorp¬ 

tion  events  at  three  magnetic  latitudes  -  II.  Differences  between  conjugate 
regions,  Planet.  Space  Sci.  1^1571. 

15.  Foppiano,  A.  (1975)  CCIR,  IWP  611,  Docs.  3  and  10. 

16.  Elkins,  T.  J.  (1972)  A  Model  of  Auroral  Substorm  Absorption,  AFCHI.-TH-72- 

0413,  AD  A749859. 
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of  the  auroral  oval.  Auroral  sub-storms  occur  sporadically  and  with  varying  in¬ 
tensity.  It  is  reasonable  to  represent  auroral  absorption  as  a  statistical  random 
variable,  recognizing  that  its  effects  are  superimposed  on  the  "normal"  absorption 

component.  When  no  auroral  substorm  is  in  progress,  the  residual  absorption 

1 2 

present  is  that  predicted  by  the  Bradley-George  model. 

In  the  present  study,  auroral  absorption  will  be  represented  by  its  median 
value,  derived  as  a  function  of  geomagnetic  latitude  and  time  from  the  references 
cited  in  the  previous  paragraph.  Since  auroral  absorption  is  ordered  in  the  cor¬ 
rected  geomagnetic  coordinate  system,  an  additional  complexity  is  introduced  in 
including  it,  since  the  normal  D- region  absorption  is  primarily  a  function  of  geo¬ 
graphic  latitude  and  time.  Figure  5  shows,  in  the  CGS,  the  percentage  of  time  that 

13 

one  way,  vertical  incidence  absorption  at  30  MHz  exceeds  one  dB.  The  contours 
are  labelled  in  percent.  Note  that  maximum  auroral  absorption  occurs  between 
08  and  09  CGT  and  at  65°  CGL.  In  order  to  use  this  data  in  computing  oblique 
incidence  HF  clutter  absorption,  median  values  for  the  statistical  distributions  of 
absorption  values  need  to  be  derived,  and  appropriate  scaling  needs  to  be  applied 
for  frequency  and  obliquity  of  propagation. 

2.4  Magnetic  Field  Model 

Several  models  exist,  with  varying  degrees  of  complexity,  to  represent  the 

structure  of  the  geomagnetic  field  in  those  dimensions.  For  an  application  such  as 

the  present  one,  in  which  it  is  necessary  to  locate  the  "orthogonality  point"  to  an 

accuracy  considerably  better  than  100  km,  the  simplest  model  (that  of  an  eccentric 

dipole)  is  much  too  inaccurate.  The  two  models  which  were  considered  for  this 

study  were  the  International  Geomagnetic  Reference  Field  -  IGRK**'  and  the 
ii  19 

" Jensen-Cain"  model.  These  are  both  higher  order  spherical  harmonic  expan¬ 
sions  of  the  geomagnetic  potential;  the  IGKF  is  an  8th  order  expansion  and  the 
Jensen-Cain  model  is  a  6th  order  expansion.  Other  expansions  have  been  proposed 
having  even  a  greater  number  of  spherical  harmonic  terms  but  the  improved  accu¬ 
racy,  if  any,  inherent  in  their  use  does  not  warrant  the  additional  computational 
complexity. 


17.  1AGA  Commission  2,  Working  Group  4  (1969)  International  geomagnetic 

reference  field,  1965.  0,  .1.  Gcophys.  lies.  £4:4407. 

18.  Mead,  G.  D.  <  1970)  International  geomagnetic  reference  field  I  963.  0  in  dipole 
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Figure  5.  Contours  of  Vertical  Auroral  Absorption  at  30  MHz  (Hartz  and 
Montbriand).  Numbers  represent  percentage  of  time  absorption  exceeds 
1  dB 


In  order  to  compare  the  IGRF  and  Jensen-Cain  models  for  the  auroral  back- 
scatter  application,  a  series  of  ray  tracings  was  performed  using  an  actial  ionos¬ 
pheric  model,  and  computing  the  geomagnetic  aspect  angle  (I//)  for  both  magnetic 
field  models.  When  the  results  were  averaged  over  a  range  interval  from  the  radar 
of  200-1100  km,  a  root-mean- square  difference  in  the  computed  values  of  for 
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both  models  was  found  to  be  0.  05  degrees.  This  value  is  sufficiently  small  so  that 
no  further  effort  was  judged  necessary  in  establishing  which  of  the  two  models  is 
the  more  accurate  for  this  application.  Because  of  its  greater  simplicity  and  for 
other  reasons  of  convenience,  the  Jensen-Cain  model  was  chosen  to  represent  the 
geomagnetic  field  in  this  study. 

As  a  matter  of  interest,  a  comparison  was  made  between  the  Jensen-Cain 
model  and  the  eccentric  dipole  model.  The  angle  \f/  was  computed  at  a  radar  lo¬ 
cation  (61°  latitude,  0°  longitude)  for  zero  elevation  take-off  angle,  as  a  function 
of  azimuth.  Figure  6  shows  the  results  of  this  comparison.  Typical  values  of 
several  degrees  for  the  discrepancy  in\p  are  apparent.  The  difference  in  computed 
magnetic  inclinations  at  the  radar  is  2.  3°  and  the  difference  in  computed  declina¬ 
tions  is  13°.  It  is,  therefore,  clear  that  the  eccentric  dipole  model  is  inadequate 
for  the  present  application. 


AZIMUTH 


Figure  6.  Comparison  Between  Jensen-Cain  Multipole  Geomagnetic  Field  Model 
and  Dipole  Model,  Showing  Angle  Between  Horizon  Ray  From  a  Site  at  <61°  N; 

0°  E)  and  Geomagnetic  Field  Vector  as  a  Function  of  Azimuth.  Magnetic 
inclination  and  declination  for  the  site  are  shown  for  each  model 
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3.  IONOSPHERIC  IRREGULARITY  MORPHOLOGY 


A  reliable  predictive  model  for  auroral  radar  clutter  must,  of  necessity,  take 
into  account  the  spatial,  temporal,  seasonal,  and  magnetic  activity  dependences  of 
ionospheric  irregularity  occurrence.  Unfortunately,  these  have  not  been  sufficient¬ 
ly  well  studied  by  the  backscatter  radar  technique  to  permit  construction  of  a  model 
based  on  data  of  the  most  relevant  kind— namely  HF  backscatter.  On  the  other 
hand,  the  study  of  ionospheric  irregularities  has  always  been  a  prominent  part  of 
ionospheric  radio  physics  and  an  extensive  data  base  exists,  accumulated  over 
many  years  by  a  variety  of  techniques.  The  most  comprehensive  component  of  this 
data  base  is  that  collected  by  the  world-wide  network  of  vertical  incidence  ionos¬ 
pheric  sounders.  Figure  3  shows  a  map  of  the  northern  hemisphere,  with  the 
auroral  oval  superimposed  and  with  the  locations  of  high  latitude  vertical  incidence 
sounding  stations  marked.  It  can  be  seen  that  this  network  has  achieved  reason¬ 
ably  well  distributed  coverage  of  the  high  latitudes.  A  similar,  but  sparser  network 
has  operated  in  the  southern  high  latitudes  and  some  of  this  data  has  been  incorpor¬ 
ated  into  the  model  as  well  as  the  northern  data. 

3.1  E-region  Irregularities 

Small  scale  irregularities  in  the  E-region  give  rise  to  the  phenomenon  of 

sporadic-E  reflections,  observed  in  routine  vertical  incidence  sounding  of  the 

ionosphere.  The  most  useful  study  of  sporadic-E  morphology  is  that  of  Besproz- 
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vannaya  and  Shchuka  who  developed  percentage  occurrence  statistics  in  the  CGS. 
Specifically,  these  authors  presented  the  percentage  of  time  that  the  maximum  fre¬ 
quency  of  high-latitude  sporadic-E  reflections  exceeded  3  MHz,  as  a  function  of 
geomagnetic  latitude,  geomagnetic  time,  season,  magnetic  activity,  and  sunspot 
number.  In  retrospect,  the  choice  of  such  a  low  maximum  frequency  as  3  MHz 
was  somewhat  unfortunate  since  a  "saturation"  effect  of  occurrence  values  close 
to  100  percent  is  often  in  evidence.  Despite  this  minor  shortcoming,  the  model 
for  E-region  auroral  clutter  is  based  on  the  above-mentioned  study.  Implicit  in 
this  procedure  is  the  assumption  that  the  same  ionospheric  irregularities  give  rise 
to  both  vertical  incidence  sporadic-E  and  oblique  incidence  HF  auroral  clutter. 

This  assumption  is  regarded  as  valid  in  view  of  the  similar  frequency  ranges  in¬ 
volved  in  observing  both  phenomena.  It  would  be  more  in  question  if  the  model 
were  intended  to  characterize  VHF  or  IJHF  auroral  clutter,  since  the  irregularity 
scale  size  regimes  in  question  would  then  be  different  by  several  orders  of 


21.  Hcsprozvannaya,  A.S. ,  and  Shchuka,  T.  I.  ( 197  1)  Geomag,  and  Aeron.  6^878; 
and  (1972)  Space-time  distribution  of  E„  formations  associated  with  visible 
auroral  forms.  Geomag.  and  Aeron,  12:401. 


magnitude;  since  quite  different  plasma  instability  processes  might  be  dominant 
at  low  and  high  wave  numbers,  the  assumption  of  such  a  direct  relationship  between 
sporadie-E  and  UHF  backscatter  clutter  is  questionable. 

The  model  for  E- region  auroral  clutter  incidence  exhibits  the  following  general 
properties: 

(1)  A  strong  tendency  for  geomagnetic  latitude  control, 

(2)  A  strong  midnight  maximum,  in  winter  especially,  which  is 
assumed  to  be  predominantly  direct-E.  This  maximum  occurs 
within  the  auroral  oval  boundaries. 

(3)  Strong  maxima  in  occurrence  at  lower  latitudes  (~60°)  at  04 
and  19  MLT.  These  maxima  are  assumed  to  correspond  to 
diffuse  radar  returns, 

(4)  A  pronounced  minimum  in  occurrence  frequency  at  65°  CGL, 

(5)  A  pronounced  dependence  of  occurrence  frequency  on  magnetic 
activity, 

(6)  A  virtual  independence  of  occurrence  frequency  on  sunspot 
number,  after  the  solar  cycle  dependence  of  geomagnetic 
activity  has  been  removed, 

(7)  A  tendency  for  low  latitude  (~60°)  occurrence  frequencies 
to  be  greater  in  summer  than  in  winter. 

Many  of  the  features  exhibited  by  the  sporadic- E-based  models  have  also  been 

observed  in  backscatter  experiments.  Among  the  studies  of  diurnal  variations  are 
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those  of  Unwin,  Presnell  et  al,  Hellgren  and  Meos,  Brooks,  Czechowsky 
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et  al,  and  Sprenger  and  Glode.  Several  of  these  studies  confirm  some  aspects 
of  the  latitudinal  distribution  although  the  requirements  for  magnetic  aspect  angle 
are  usually  too  restrictive  to  permit  reliable  extraction  of  latitudinal  dependence 
in  a  UHF  radar. 

3.2  F-region  Irregularities 

Small  scale  irregularity  of  the  F-region  polar  ionosphere  is  a  very  common 

phenomenon— in  fact,  considering  that  the  F-region  may  be  defined  as  extending 

from  an  altitude  of  ~  170  km  to  one  of  -  1000  km,  small  scale  irregularities  in  this 

region  may  be  regarded  as  a  permanent  feature  of  the  high-latitude  ionosphere. 
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Aarons,  using  the  satellite  scintillation  technique,  reports  that  the  F-region 
ionosphere  across  the  polar  cap  (within  the  equatorial  boundary  of  the  auroral  oval, 
approximately)  is  permanently  filled  with  irregularities.  From  the  viewpoint  of 
HF  backscatter,  this  observation  is  not  as  relevant  as  for  satellite  communications, 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  87.  ) 
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for  example,  since  only  irregularities  on  the  underside  of  the  F-layer,  below  the 
F-peak,  would  normally  affect  an  HF  system.  The  model  adopted,  in  this  report, 
for  F-region  irregularity  will  be  based  on  observations  of  the  bottomside  ionos¬ 
phere,  by  vertical  incidence  sounders.  A  HF  radar  sees  F-region  irregularities 
as  "slant-F"  clutter,  named  for  its  characteristic  appearance  on  oblique-incidence 
ionograms. 

Slant-F  clutter-producing  irregularities  can  be  assumed  to  fill  the  radar  beam 

completely;  these  are  the  same  irregularities  which  produce  spread-F  on  vertical 
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incidence  ionosondes.  That  this  is  so  was  demonstrated  by  Bates,  among  others. 
If  the  F-region  irregularities  in  question  extend  over  the  radar  location  and  down 
range,  slant-F  is  seen  continuously  in  range  at  all  frequencies  which  do  not  pene¬ 
trate  the  ionosphere.  The  intensity  of  slant-F  clutter  is  proportional  to  the  illumina 

ted  volume  (V  )  within  the  irregularity  region,  and  to  the  volumetric  cross  section 
s 

of  the  irregularities.  At  high  geomagnetic  latitudes,  where  the  magnetic  field  lines 
are  nearly  vertical,  a  simple  rule  of  thumb  is  that  slant-F  clutter  occurs  at  about 
half  the  range  of  ground  clutter  as  seen  by  the  sky-wave  (for  the  same  take-off 
angle). 

3.  2.  1  MORPHOLOGY  OF  SLANT-F 

Since  slant-F  and  spread-F  are  produced  by  the  same  ionospheric  irregulari¬ 
ties,  it  is  possible  to  deduce  the  spatial  and  temporal  properties  of  slant-F  from 
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those  of  spread-F.  Spread-F  has  been  extensively  studied  on  a  global  basis. 

In  general,  the  following  properties  may  be  assigned  to  high  latitude  spread-F: 

(1)  Latitudinal  dependence:  Controlled  principally  by  geomagnetic 
latitude,  with  the  probability  of  spread-F  occurrence  increasing 
towards  high  latitudes, 

(2)  Temporal  dependence:  A  maximum  in  probability  of  occurrence 
is  found  in  the  auroral  F-region  at  around  02-04  magnetic  local 
time  (MLT).  A  second  maximum  was  reported  by  several 
authors  to  be  permanently  located  over  Northeastern  Canada, 
but  it  is  probable  that  this  is  associated  with  the  day-side 
auroral  oval  zone  of  low  energy  particle  precipitation  (see 
Appendix  A), 
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(3)  Seasonal  dependence:  Maximum  probability  of  spread-F  occurrence 
occurs  in  winter,  with  the  minimum  in  summer.  This  is  most 
likely  a  solar  illumination  phenomenon,  since  it  seems  that  high- 
latitude  spread-F  is  sharply  reduced  by  the  stabilizing  influence 

of  ionizing  solar  radiation.  During  the  long  winter  night,  spread-F 
is  a  virtually  permanent  feature  of  the  high-latitude  ionosphere.  In 
day-time,  spread-F  is  normally  only  seen  in  the  vicinity  of  the  day- 
side  auroral  zone  (projection  of  the  magnetospheric  cleft), 

(4)  Magnetic  activity  dependence;  At  high  latitudes,  there  is  a  pronounced 

tendency  for  the  occurrence  of  spread-F  to  increase  with  increasing 
30 

magnetic  activity.  This  tendency  is  more  apparent  at  the 

equinoxes  than  at  the  solstices,  when  a  less  pronounced  correlation 

with  magnetic  activity  is  observed.  The  equatorward  boundary  of 

the  spread-F  region  moves  to  lower  latitudes  with  increasing 
34 

geomagnetic  activity. 

3.  2.  2  OBSERVATIONS  OF  SLANT-F 

or  op  0*7 

The  slant-F  phenomenon  has  been  studied  by  a  number  of  authors.  '  * 

All  of  these  authors  note  the  solar  zenith  angle  control  of  the  occurrence  frequency 
of  slant-F.  Baggaley,  working  at  Sheffield,  England,  has  presented  the  most 
thorough  analysis  of  this  aspect  of  slant-F.  He  showed  that  the  onset  of  slant-F 
is  most  likely  to  occur  when  the  solar  zenith  angle  is  between  90-110°.  From 
Figures  11  and  12  of  Peterson  et  al  a  similar  conclusion  may  be  drawn,  assum¬ 
ing  that  the  echo  regions  observed  by  them  were  about  1000  km  to  the  north.  Bagga¬ 
ley  notes  that  the  mean  duration  of  slant-F  events  is  1.  5  hr,  whereas  Peterson  et 
al  indicate  that  many  events  persist  throughout  the  night,  until  sunrise  in  the  ionos¬ 
phere.  The  difference  between  these  two  observations  is  that  Baggaley's  were  made 
at  sunspot  minimum  at  a  northerly  latitude,  while  those  of  Peterson  et  al  were  made 
at  a  higher  sunspot  number  and  at  a  more  southerly  latitude.  After  sunset,  the 
17  MHz  transmitted  beam  of  Baggaley  penetrated  the  ionosphere  as  the  latter's 
critical  frequency  diminished,  whereas  the  ionosphe re  to  the  north  of  Stanford, 
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where  the  observations  of  Peterson  et  al  were  made,  was  sufficiently  strong  to 
reflect  their  17.  3  MHz  transmissions  throughout  most  nights. 

All  three  authors  mentioned  above,  note  that  the  F- region  slant  echoes  display 
a  range  dependence  exhibiting  an  abrupt  onset  at  a  short  range  which  may  vary 
considerably.  This  observation  is  also  related  to  ionospheric  penetration,  with  the 
abrupt  inner  boundary  of  the  echo  region  corresponding  to  one  half  of  the  classical 
"skip  distance"  for  ground  to  ground  propagation.  All  three  authors  note  the  aspect 
sensitivity  of  the  echoes,  but  the  observations  of  Baggaley  are  especially  interest¬ 
ing  because  they  permit  a  quantitative  estimate  to  be  made  of  this  aspect  sensitivity. 
Baggaley.  in  his  Figure  7,  gives  the  dependence  of  echo  activity  on  magnetic  aspect 
angle,  based  on  the  azimuthal  variation  of  the  minimum  echo  range.  He  shows  the 
echo  strength  dropping  below  noise  level  (presumably)  at  10  degrees  off-orthogonal¬ 
ity,  and  the  maximum  number  of  echoes  at  5  degrees.  By  relating  these  results  to 

the  probability  distribution  of  F- region  echoes  recorded  in  the  Polar  Fox  II  experi- 
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ment,  taking  into  account  the  different  sensitivities  of  the  two  equipments,  a 
figure  of  about  10  dB/degree  off  perpendicularity  is  reached. 

Baggaley  shows  that  radial  velocity  of  F- region  scatterers  may  lie  between 
0  and  150  msec,  with  a  mean  value  of  65  m/sec.  This  observation  made  in  the 
United  Kingdom,  is  relevant  for  a  magnetic  latitude  of  about  63°,  at  night. 

3.3  Plasma  Motion 

The  Doppler  shift  experienced  by  a  radar  signal  backscattered  from  ionos¬ 
pheric  irregularities  is  a  result  of  the  motion  of  the  plasma  in  which  these  irregu¬ 
larities  arc  embedded,  together  with  any  component  of  motion  of  the  irregularities 
relative  to  the  plasma.  The  motion  of  ionospheric  plasma  in  the  earth's  polar 
regions  is  predominantly  due  to  the  “lectrodynamic  effects  induced  by  the  mag- 
netospheric  electric  potential  di~4ribution.  This  takes  the  form  of  a  linear  gradient 
of  electric  potential  in  the  east-we. J  direction  across  the  earth's  magnetospheric 
tail  (see  Section  1.  3).  Figure  2  illustrates  the  convection  pattern  induced  by  this 
electric  field  in  the  equatorial  plane  of  the  magnetosphere.  The  convection  pattern 

in  the  ionosphere  may  be  visualized  by  mapping  the  magnetospheric  convection 
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pattern  down  into  the  ionosphere  along  the  geomagnetic  field  lines.  Maeda  has 
constructed  a  model  for  the  E-  and  F-region  plasma  drifts,  based  on  a  realistic 
representation  of  the  cross-tail  electric  field  and  the  electrical  conductivity  of  the 
ionosphere.  Maeda's  numerical  estimates  of  ionospheric  plasma  velocities  are 
in  general  agreement  with  many  measured  values. 
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The  basic  pattern  of  ionospheric  plasma  drift  is  a  "two-cell"  pattern  in  which 
counter-rotating  flows  around  the  periphery  of  the  auroral  oval  close  in  an  anti- 
sunward  direction  across  the  polar  cap.  Figure  7  contains  an  additional  plasma 
drift  component  due  to  co- rotation  of  the  earth's  atmosphere  with  the  solid  earth 
(see  Section  1.  3)  and  is  thus  the  pattern  which  would  be  observed  from  a  non¬ 
rotating  frame  of  reference,  for  example,  a  satellite.  This  co-rotational  compo¬ 
nent  is  relatively  small  compared  to  the  cross-tail  magnetospheric  electric  field 

component.  This  pattern  has  been  substantiated  and  confirmed  by  a  variety  of 

4 1 

satellite  and  ground-based  measurements.  Heelis  et  al  using  the  Atmospheric 

Explorer-C  satellite  at  an  altitude  of  250  km,  measured  plasma  velocities  in  the 

region  of  the  dayside  cleft  and  reported  velocities  in  excess  of  l  km /sec,  exhibiting 

a  reversal  near  the  noon  meridian,  which  is  characteristic  of  the  two-cell  pattern. 
42 

Tsunoda  et  al  used  data  from  a  ground-based  backscatter  sounder  in  the  polar 

cap  and  confirmed  the  two-cell  pattern  by  noting  the  azimuthal  pattern  of  slant-E 

echoes.  Similar  inferences  can  be  drawn  from  the  studies  of  polar  cap  sporadic-E 
43 

by  Taieb  and  from  measurements  of  polar  cap  auroral  backscatter  made  during 

the  Polar  Cap  III  experiment.  A  useful  review  of  measurements  of  plasma  drift 

44 

velocities  at  high  latitudes  is  contained  in  the  works  of  Doupnik  et  al,  and 
45 

Doupnik. 
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Figure  7.  A  Model  of 
F-layer  Plasma  Drift, 
Including  Effect  of  Earth 
Rotation  (after  Heppner, 
1977).  The  drift  speed  is 
proportional  to  the  density 
of  flow  lines,  being  about 
0.  5  km /sec  near  the  pole. 
The  coordinate  system  is 
..  the  CGS.  The  auroral 
oval  (Q  =  3)  is  shown 
dotted  and  the  region  of 
magnetiospheric  cleft 
precipitation  near  noon 
is  shown  dashed 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  87.) 
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The  velocity  pattern  for  E-region  plasma  differs  significantly  from  the  E  -  region 
pattern.  Due  to  the  very  low  E-region  collision  rates,  the  llall  conductivity  of  the 
K-region  is  negligibly  small,  whereas  in  the  E-region  it  is  very  high.  Thus  the 
motion  of  plasma  in  the  F- region  is  in  the  direction  of  the  imposed  electric  field, 
whereas  in  the  E-region,  there  is  a  large  component  orthogonal  to  the  electric 
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field.  The  analysis  of  Maeda  has  taken  into  account  these  and  other  factors  and 
will  be  used  as  the  basis  of  the  model  to  be  developed  for  radar  Doppler  shifts 
in  Section  5. 

3.4  Summary  of  Auroral  Kadar  Experimrnla 

A  wealth  of  experimental  evidence  has  been  collected  over  a  span  of  three 

decades  on  the  subject  of  the  radio  aurora,  in  spite  of  this  large  data  base,  many 
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aspects  of  this  complex  phenomenon  remain  obscure.  Unwin  published  a  general 
survey  of  the  morphology  of  radar  aurora  including  its  diurnal  and  seasonal  varia¬ 
tions,  noting  that  E- region  radio  aurora  manifests  itself  in  two  distinctly  different 
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forms  — diffuse  and  discrete.  Unwin  and  llaggaley  have  published  an  updated 
review  of  the  E-region  radio  aurora,  and  a  great  many  authors  have  reported 
experimental  studies  performed  with  a  variety  of  instruments.  An  exhaustive 
survey  of  this  work  would  be  far  too  lengthy  to  be  included  in  this  section,  and  only 
major  findings  will  be  discussed.  In  particular,  since  primary  interest  in  this  work 
is  in  connection  with  HE  radio  aurora,  special  attention  will  be  given  to  observa¬ 
tions  made  at  HE  and  low  VHE. 

Diffuse  aurora  exhibits  little  structure  in  both  spatial  and  Doppler  domains, 
and  is  found  to  exhibit  a  close  relationship  to  the  eastward  auroral  electrojet, 
occurring  equatorward  of  the  statistical  auroral  oval  in  the  late  afternoon  and  even¬ 
ing.  Discrete  radio  aurora  is  found  predominantly  in  the  hours  around  midnight  and 
shortly  after,  and  exhibits  sharply  defined  structure  in  space  and  velocity.  Dis¬ 
crete  aurora  is  associated  with  the  so-called  "break-up"  phase  of  the  visual  aurora. 
While  discrete  auroral  returns  are  rare  in  the  late  afternoon  and  evening,  diffuse 

returns  are  observed  at  all  times,  though  in  combination  with  discrete  returns  after 
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midnight.  Greenwald  et  al  have  demonstrated  convincingly  that  the  diffuse 

auroral  backscatter  is  collocated  with  auroral  currents.  The  observations  of 
46 

Unwin  demonstrate  that  the  currents  associated  with  the  aurora  at  CGI.  ~60°  are 


46.  Unwin,  R.  S.  (1966)  The  morphology  of  VHE  radio  aurora  at  sunspot  maxi¬ 

mum  II.  The  behavior  of  different  echo  types,  J.  Atmos.  Terr.  I’hys. 
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frequently  intense  enough  to  produce  strong  radio  aurora  at  that  latitude.  Many 
49  50  51 

other  experiments,  '  ’  also  have  demonstrated  appreciable  radar  backscatter 

from  sub-auroral  latitudes. 

The  diurnal  variation  of  auroral  backscatter  exhibits  two  maxima,  correspond¬ 
ing  to  the  diffuse  and  discrete  echo  types.  Summarizing  the  results  of  many  experi¬ 
ments,  too  numerous  to  cite  individually,  the  evening  maximum  of  diffuse  aurora 
occurs  between  16  and  19  CGT,  and  the  "midnight"  maximum  of  discrete  echoes 
occurs  between  23  CGT  and  02  CGT.  The  intensity  of  auroral  backscatter  also 
exhibits  periodicities  on  a  long  time  scale  (seasonal  and  solar  cycle).  The  seasonal 
variation  exhibits  maxima  at  the  equinoxes,  with  a  minimum  in  summer  and  a  still 
deeper  minimum  in  winter.  In  terms  of  probability  of  occurrence,  the  equinoctial 
maxima  are  roughly  twice  as  great  as  the  winter  minimum.  The  solar  cycle  de¬ 
pendence  shows  a  positive  correlation  between  radio  auroral  activity  and  sunspot 
number,  with  the  former  lagging  the  latter  by  1-2  years.  These  seasonal  and  solar 
cycle  periodicities  are  consistent  with  the  generation  of  auroral  irregularities  by 
ionospheric  currents,  being  very  similar  to  the  corresponding  long  term  periodici¬ 
ties  in  magnetic  activity  index  (see  Section  6).  In  fact,  as  many  authors  have 

demonstrated,  auroral  radar  intensity  is  strongly  correlated  with  geomagnetic 
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activity.  Unger  et  al  obtained  better  correlation  between  auroral  radar  activity 
and  a  composite  "solar-planetary"  index  (the  product  of  K^  and  sunspot  number) 
than  with  either  index  separately. 

The  observations  noted  thus  far  were  made  at  VHF  and  UHF  and  relate  to  the 

E-region,  since  geometric  constraints  make  F-region  irregularities  inaccessible 

to  line-of-sight  radars.  In  order  to  observe  ionospheric  backscatter  from  the  F- 

region,  observations  at  HF  are  necessary,  utilizing  ionospheric  refraction  to 

achieve  magnetic  orthogonality  at  the  required  altitudes.  The  most  comprehensive 
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HF  auroral  radar  experiment  yet  performed  was  the  Polar  Fox  II  experiment  in 
Northern  Maine.  The  results  of  this  experiment  will  be  used  extensively  in  develop¬ 
ing  the  clutter  model  in  this  report,  F-region  backscatter  was  found  by  Polar  Fox  II 
to  be  virtually  always  present,  during  nighttime  at  some  range.  Another  thorough 

49.  Flood,  W.A.  (1960)  VHF  auroral  backscatter  measurements,  J.  Geophys.  Res. 
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50.  Gartlein,  C.  W. ,  Sprague,  G. ,  and  Waag,  R.  C.  (1960)  41  Mc/s  IGY  auroral 

radar  at  Ithaca,  New  York,  J.  Geophys.  Res.  65: 2255. 

51.  Watkins,  C.D.  (196 1)  Auroral  radio-echoes  and  magnetic  disturbance, 

J.  Atmos.  Terr.  Phys.  20:131. 
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study  of  HF  auroral  backscatter  was  conducted  at  Ava,  New  York,  at  CGL  =  56*. 
The  maximum  occurrence  of  auroral  associated  F-region  backscatter  was  found  to 
occur  between  20  MLT  and  08  MLT  in  both  winter  and  summer.  The  probability 
of  F- region  backscatter,  at  03  MLT  at  the  ionospheric  intersection  point  (~67*  CGL) 
was  over  90  percent.  A  high  power  backscatter  experiment  conducted  in  the  south 
of  England  produced  similar  results.  The  probability  of  F-region  backscatter  in 
this  experiment  was  over  80  percent,  and  the  amplitude  of  the  returns  from  the 
auroral  region  exceeded  that  from  the  temperate  latitude  ionosphere  by  30-50  dB. 
The  F-region  backscatter  was  found  to  have  a  sharp  leading  edge  and  a  gradual 
trailing  edge,  and  the  apparent  scatter  cross  section,  in  the  range  8-10  MHz  was 
found  to  increase  with  decreasing  frequency.  In  one  particular  instance,  the  sharp 
leading  edge  of  F-backscatter  was  observed  at  7  00  km,  when  the  operating  fre¬ 
quency  was  10.  1  MHz  and  the  F-region  critical  frequency  was  3  MHz.  An  analysis 
performed  in  Section  4.  6  will  show  that  this  observation  is  consistent  with  refracted 
propagation  through  an  irregular  F-region.  The  amplitude  of  F-region  clutter  was 
found  to  be  strongly  correlated  with  magnetic  activity. 
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Other  experiments  performed  at  HF  include  those  of  Hower  et  al  ,  Weaver, 

56  25  35  36  57 

Swenson,  Brooks.  Baggaley,  Peterson  et  al,  Malik  and  Aarons, 
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Czeckowsky  et  al,  and  Sprenger  and  Glode.  In  general  these  experiments  all 

show  F-region  auroral  related  backscatter  to  be  a  highly  common  phenomenon, 

especially  at  night  and  increasing  in  intensity  with  increasing  geomagnetic  activity. 

Another  important  fact  which  these  experiments  revealed  is  that  the  characteristics 

of  the  ionospheric  propagation  path  between  radar  and  irregularities  are  a  factor 

of  major  importance  in  determining  the  probability  of  F-region  backscatter. 

58 

Recently,  Oksman  et  al  have  conducted  a  comparison  between  VHF  auroral 
measurements  and  HF  backscatter  observations.  They  have  shown  that  the  back¬ 
scatter  cross  section  of  the  radio  aurora  follows  a  wavelength  cubed  low  in  this 
frequency  range.  These  workers  found  that  F-region  radio  aurora  often  occurred 
on  the  same  magnetic  L-shells  as  E- region  radio  aurora  and  often  occurred  in  two 
different  geomagnetic  latitude  regions— one  associated  with  the  auroral  oval  and  a 

54.  Katz,  A.H.  (1971)  HF  Auroral  Backscatter  Study,  AVCO  Kept.  AVSD-0378- 
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radar  echoes  and  high  latitude  spread- F  measurements,  J.  Geophys.  Res. 
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second  with  the  poleward  wall  of  the  ionospheric  trough.  This  second  group  of 
trough-related  echoes  appears  only  at  sunspot  minimum  and  is  absent  at  sunspot 
maximum. 

Many  workers  have  studied  the  Doppler  characteristics  of  the  radio  aurora  over 
many  years.  A  controversy  exists  over  whether  or  not  auroral  plasma  irregulari¬ 
ties  travel  at  the  local  ion-acoustic  velocity,  or  at  the  electron  drift  velocity  or  at 
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some  combination  of  these  two.  Measurements  at  UHF  tend  to  support  the  first 

of  these  alternatives,  but  measurements  at  lower  frequencies  (for  example,  at  50 

60  6 1 
MHz)  tend  to  support  the  second  alternative.  A  recent  paper  by  Greenwald  et  al 

has  discussed  this  "velocity  question".  In  developing  the  velocity /Doppler  model 

in  this  report,  it  will  be  assumed  that  the  apparent  velocity  of  auroral  scatterers, 

in  the  HF  range,  is  that  of  the  electron  drift  velocity. 

3.5  Altitude  of  the  Radio  Aurora 

Many  studies  have  established  that  the  E- region  radio  aurora  is  confined  to  a 

22 

narrow  range  of  altitudes,  centered  at  about  110  km.  For  example,  Unwin 

deduced  a  height  of  110.  9  ±  3.  1  km,  with  a  vertical  thickness  of  the  radio  scatter- 
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ing  layer  of  ~5  km.  Abel  and  Newell  reported  an  altitude  of  ~  1 10  km  with  thick- 
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ness  of  10-15  km.  Barber  et  al  reported  an  altitude  of  110  ±  1  km  with  a  thickness 

52 

of  a  few  kilometers.  Unger  et  al  report  an  altitude  of  110-120  km  with  a  thickness 
of  6-8  km.  The  altitude  of  F-region  radio  scatterers  in  the  auroral  region  is,  on  the 
other  hand,  much  less  well-defined.  Irregularities  almost  certainly  exist  through¬ 
out  the  entire  altitude  range  of  the  F-region,  from  about  180-200  km  and  upwards. 
The  intensity  of  these  F-region  irregularities  is  roughly  proportional  to  the  ambient 
electron  density,  so  that  the  most  intense  irregularities  may  be  expected  near  the 
peak  of  the  F-layer  (typically  250-350  km). 

In  this  report,  the  altitude  of  the  E-region  radio  aurora  will  be  taken  as  110  km. 
The  average  vertical  thickness  of  the  scattering  layer  will  be  taken  as  ±  5  km  be¬ 
tween  -10  dB  points,  that  is,  between  altitudes  at  which  the  backscatter  volumetric 
cross  section  is  10  dB  lower  than  its  value  at  110  km.  Where  it  is  necessary  to 

59.  Abel,  W.  J. ,  and  Newell,  R.  E.  (1969)  Measurements  of  the  afternoon  radio 
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assume  an  altitude  for  F- region  auroral  scatterers,  this  will  be  taken  as  any  alti¬ 
tude  above  180  km.  At  night,  the  peak  F- region  volumetric  backscatter  cross 
section  will  be  assumed  to  occur  at  300  km. 


4.  AURORAL  CLUTTER  MODEL 

As  a  first  step  towards  a  model  predicting  the  magnitude  of  auroral  backscatter. 
the  morphological  data  discussed  in  Section  3  has  been  used  to  construct  a  model 
for  the  distribution,  in  space  and  time,  of  the  probability  of  irregularity  occurrence. 
This  model  will  subsequently  be  translated  into  a  model  for  the  auroral  backscatter 
magnitude. 

4.1  E-  Region  Probability  of  Occurrence 

The  space-time  distribution  of  the  probability  that  fE  exceeds  3  MHz,  deduced 
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by  Besprozvannaya  and  Shchuka,  referred  to  in  Section  3.  1,  is  used  as  the  basis 
for  the  quantitative  model  for  E-region  irregularity  occurrence.  These  authors 
reported  that  auroral  sporadic-E  appeared  to  be  confined  to  two  ring-shaped  maxima 
encircling  the  geomagnetic  pole,  separated  by  a  minimum  at  about  65°  geomagnetic 
latitude.  If  the  probability  that  fEg  -  3  MHz  is  represented  by  the  product  of  a  func¬ 
tion  of  geomagnetic  latitude  and  one  of  geomagnetic  time,  it  may  be  written  as 


Pj(E)  =  Pl(4)  •  Pj(T) 


for  the  high  latitude  ring,  and 


P2<E)  =  p2<<w  •  P2<T) 


for  the  low  latitude  ring  where 


p  j  (<*»)  =  exp 


(*-  A) 


A  =  72-  1.5  (K  -2)  +  2  cos  -pr  (T  -  18) 

p  1 2 


p  j(T)  =  0.  5  +  0.  3  cos  (T  -  1) 


p2(*)  =  exp  - 


<*  -  63)* 


and,  letting  G  =  min(D;  365  -  D) 


P2(T)  =  log  (Kp  +  l)2  [0,  50  -  0.  125  cos  yj  (T  -  3)  •  (1  -  cos 

-  0.  5  exp  (-0.  1(T  -  1 9) 2  cos 

-  0.45  exp  -  (6.7  X  10'5G2)  exp  -  (0.25  (T  -  10>2)] 

and  Kp  is  planetary  magnetic  index.  There  is  also  a  small  peak  in  sporadic-E 
occurrence  at  the  geomagnetic  pole,  which  can  be  represented  by 

P3(E)  =  0.  3  exp  -  . 

In  these  expressions,  the  units  for  <t  are  degrees  and  those  for  T  are  hours.  The 
Gaussian  functions  are  chosen  for  their  convenient  asymptotic  properties,  and 
have  no  other  special  significance.  The  probability  that  fE  .  -  3  MHz  is  represented 
by  the  sum  of  the  three  components  Pj(E)  +  p2<E)  +  Pj(E). 

The  magnetic  activity  dependence  of  the  low  latitude  ring,  contained  in  p.,(T), 
is  an  attempt  to  fit  approximately  the  published  data  by  an  arbitrary  mathematical 
function  for  purposes  of  convenient  numerical  computation. 

A  computer  program  has  been  written  to  compute  P(E>  for  arbitrary  location 
and  time  and  magnetic  activity.  Figure  8  (a - f )  shows  the  distribution  in  C'GS  of 
the  resultant  probability  of  E- region  irregularity  for  summer,  winter  and  equinox 
seasons  and  for  two  levels  of  geomagnetic  activity  (Kp  =  2,4).  In  order  to  illustrate 
the  range  — azimuth  dependence  of  p(E)  from  a  specific  site,  f  igure  9  has  been 
produced  for  a  location  at  <t>  =  61°.  Contours  are  shown  for  summer  at  intervals 
of  0.  2  (that  is,  20  percent  probability)  for  Universal  midnight  and  K  =  2.  The 
rough  outline  of  the  auroral  oval  is  evident,  as  is  the  minimum  in  E^  occurrence 
at  65°  magnetic  latitude  and  the  peak  near  the  geomagnetic  pole. 


(a) 

Figure  8.  Contours  of  Probability  of  E-region  Irregularity  in  CG.S  at  12  UT 
for  Summer,  Equinox  and  Winter  arid  for  K  =2  and  4.  Contours  are  in 
intervals  of  0.  2  P 
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Figure  !>.  Probability  of  K- region  irregularity  as  a  Function  of  Range  and 
Azimuth  for  Site  at  61°  CGL;  June.  K  =  2,  UT  =  00.  Contours  labelled  in 
percentage  ^ 


■1.2  F-rpgion  Proliahilily  of  Occurrence 

Section  3.  2  contained  a  discussion  of  the  morphology  of  F- region  small  scale 
irregularities  in  the  polar  region  and  their  general  relationship  to  1IF  backscattor 
(slant-F).  A  model  has  been  constructed  based  primarily  on  the  references  cited 
in  that  section,  and  this  model  has  been  compared  to  available  backscattor  data  in 
the  auroral  and  sub-auroral  regions. 

The  F-region  irregularity  model  is  considered  in  two  separate  parts— the 
nighttime  auroral  maximum  and  the  day-side  magnetospheric  cleft  region.  The 
nighttime  component  expresses  the  probability  of  F-irregularity  as 

PjlF)  =  P|(4>)  •  Pj(T) 
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whe  re 


PjCtO  is  a  geomagnetic  latitude  term 
Pj(T)  is  a  geomagnetic  time  term 

and 

p,  (<*>)  =  exp  -  a[4>  -  (70.  5  -  1.  5  K  )]2 
1  P 

Pj(T)  =  exp  -  b 
whe  re 

a  =  +  0.  01178 

b  =  +  0.4834 
Kp  =  planetary  magnetic  index 

and 

<t>  is1  measured  in  radians. 

The  dayside  cleft  region  is  modelled  similarly,  as  the  product  of  a  geomagnetic 
latitude  function  and  a  geomagnetic  time  function 

P2<f'')  =  P2(<1>)  ‘  P2m 
where 

P2(<W  =  exp  -  [0.  1337<4>  -  A)2] 

P2(T)  =  exp  -  (0.  1003(T  -  12)2] 

and 

A  =  78  +  cos  <UT  -  18)  +  cos  (D  -  182.  5)  -  1.5  cos  |  (T  -  6)  . 
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(T  -  2) 

V1 


UT  =  Universal  Time 


and 


D  =  Julian  day  number. 

This  model  contains  latitudinal,  seasonal,  magnetic,  and  Universal  Time  dependen¬ 
cies.  It  is  based  on  the  work  of  Whalen  and  Pike*'^  who  measured  the  position  of 

the  equatorward  boundary  of  the  irregularity  zone  with  an  airborne  ionosonde,  and 

( j  4  t)  3 

modified  by  the  work  of  Burch  and  Oguti  and  Marubashi.  ‘  These  latter  authors 
showed,  respectively,  the  dependence  on  the  angle  between  the  earth's  magnetic- 
axis  and  the  ecliptic  plane  and  geomagnetic  latitude  and  time  dependences  of  the 
dayside  cleft  region. 

Because  of  the  pronounced  seasonal  variations  in  the  lengths  of  the  solar- 
illuminated  day  at  high  latitudes,  and  because  many  authors  (for  example,  some 
cited  in  Section  3.  2)  note  that  solar  illumination  inhibits  formation  of  F- region 
irregularities,  it  was  necessary  to  add  a  third  component  of  probability  of  F-irregu- 
larity  formation,  dependent  on  solar  zenith  angle  <\). 

Thus 

P3(F)=  Aexp-  101’ 

whe  re 

A  =  0.  2  cos  ^ ^  +  0.  7 
P3U-’)=0;  X<  90° 


and 


P3<F)  =  0.  909  (x  -  90)  •  p3(F) 


X=  101° 


for  90°  <  x  <  101°. 
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65.  Oguti,  T. ,  and  Marubashi,  K.  (1966)  Enhanced  ionization  in  the  ionospheric 

Fa  region  around  geomagnetic  noon  in  high  latitudes,  Repts.  lonos.  Spa. 
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In  these  equations  there  exist  three  domains  of  When  the  sun  does  not  illuminate 
the  lower  F- region  (x  >  101°)  we  have  a  geomagnetic  latitude  dependence  and  a 
seasonal  dependence  of  irregularity  probability.  When  the  sun  illuminates  the 
ground  (x  <  90°)  we  have  total  inhibition  of  irregularity  formation.  Between  these 
two  domains  there  is  a  linear  matching  function  ensuring  a  non-abrupt  transition. 

The  final  expression  for  the  probability  of  F-region  irregularity  formation  is 

p(F)  =  p  j(F)  +  P2<F)  +  Pg(F) 

subject  to  the  condition  that  p(F)  cannot  exceed  unity.  A  computer  program  has 
been  written  to  compute  p(F)  for  arbitrary  location  and  geophysical  variables.  Fig¬ 
ure  10  (a-f)  shows  the  distribution  in  CGS  of  the  resultant  probability  of  F- region 
irregularity  for  summer,  winter,  and  equinox  seasons  and  for  two  levels  of  geo¬ 
magnetic  activity  (K^  =  2,  1).  By  taking  cuts  through  one  of  these  contour  plots. 
Figures  11  and  12  have  been  produced,  showing  the  probability  of  irregularity  as  a 
function  of  CGL  for  00,  03,  06,  and  12  MLT,  in  summer  for  moderate  magnetic 
activity. 

In  order  to  illustrate  the  range-azimuth  dependence  of  p(F)  from  a  specific  site. 
Figure  13  has  been  produced  for  a  site  at  4  =  61°.  Contours  are  shown  for  summer, 
at  intervals  of  0.  2  (that  is,  20  percent  probability)  for  Universal  midnight.  The 
general  outline  of  the  auroral  oval  is  apparent  in  this  figure,  with  the  nighttime 
maximum  occurring  close  to  the  observing  location  and  the  daytime  C'cleft"  region) 
maximum  occurring  at  a  range  of  4000  km. 

1.3  Scatter  Volume  and  Volume  Keflertivily 

The  estimation  of  the  backscuttered  power  from  the  aurora  is  a  complicated 
procedure  in  the  case  of  an  IIF  radar.  It  is  necessary  first  to  compute  the  volume 
of  auroral  scatterers  illuminated  by  the  radar  and  the  magnetic  aspect  angle;  a 
knowledge  of  the  volume  reflectivity  of  the  auroral  scatterers,  as  a  function  of 
aspect  angle  and  frequency,  is  then  required  to  complete  the  computation.  The 
calculation  of  illuminated  volume  in  the  presence  of  substantial  ionospheric  refrac¬ 
tion  is  very  different  from  the  corresponding  calculation  at  VHF  and  I  MF  where 
straight-line  radio  propagation  can  be  assumed  with  accuracy.  Appendix  contains 
a  discussion  of  this  subject,  and  an  expression  is  derived  for  the  scatter  volume  in 
terms  of  radar  parameters  (antenna  bearmvidth,  equivalent  pulse  length)  and 
properties  of  the  l(F  ray  trajectories,  which  themselves  depend  on  ionospheric 
structure  and  radar  frequency. 


(b) 

Figure  10.  (Cont) 


45 


FUNCTION 


PROS  <SF) 


AZIMUTH 


Figure  13.  Same  as  Figure  9,  Except  for  F- region 


The  scatter  volume  for  an  HF  radar,  under  certain  conditions,  may  be  com¬ 
puted  relatively  simply  in  a  manner  analagous  to  the  procedure  for  VHF  and  U11F 
radars.  At  night  and/or  at  relatively  high  frequencies  in  the  HF  band,  the  refrac¬ 
tion  in  the  ionosphere  may  be  small.  Certainly,  under  these  conditions  there  is 
very  little  ionization  in  the  lower  ionosphere,  so  that  it  is  legitimate  to  compute 
E-region  clutter  amplitude  on  the  basis  of  straight-line  propagation.  In  this  case 
the  scatter  volume  is  determined  solely  by  radar  pulse  length  and  antenna  azimuthal 
beamwidth  (see  Figure  14). 

The  determination  of  volumetric  cross  section,  that  is,  the  effective  point 
target  cross  section  per  unit  volume  of  scatterers,  is  complicated  by  the  same 
factors  which  make  the  computation  of  scatter  volume  difficult.  Auroral  back- 
scatter  measurements  at  HF  have  never  been  translated  into  reliable  volumetric 
cross-section  measurements  because  of  the  difficulty  in  computing  the  scatter 


volume.  One  approach  to  overcoming  this  problem  is  to  extrapolate  HF  cross 
sections  from  VHF  and  U HF  measurements,  but  the  frequency  scaling  law  is  not 
known  with  reasonable  precision,  partly  due  to  the  uncertainty  whether  the  relevant 
plasma  irregularities  in  each  case  result  from  a  common  physical  process.  The 
approach  adopted  in  this  report  has  been  to  extrapolate  data  from  a  HF  radar 
experiment  to  any  location  and  time  by  use  of  the  occurrence  probability  statistics 
developed  in  this  section. 


AURORAL 

LAYER 


Figure  14.  Illustrating  Intersection  of  a  Radar  Pulse 
With  a  Thin  Auroral  Layer.  The  scattering  volume  is 
determined  by  azimuthal  beamwidth,  pulse  length,  and 
layer  thickness 


Data  from  the  Polar  Fox  II  experiment  have  been  analyzed  at  RADC'66  to  pro¬ 
vide  a  "pseudo  volumetric  cross  section."  Data  for  the  nighttime  hours  of  October 
1971  were  reduced  to  equivalent  point  target  cross  sections  by  solving  the  radar 
equation  using  Polar  Fox  II  radar  parameters.  A  representative  value  for  illumina¬ 
ted  volume  was  used,  based  on  the  geometry  of  ionospheric  ray  paths  in  an  average 
nighttime  ionosphere  (Figure  15,  Sales),  to  reduce  this  to  a  volumetric  cross 


66.  Ring,  W.,  and  Richards,  E.  (1978)  Fall  URSI  Symposium,  Boulder,  Colorado. 

67.  Sales,  G.S.  (1977)  URSI  Spring  Meeting,  Palo  Alto,  California. 
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section.  Ionospheric  absorption  was  taken  into  account  in  deducing  these  cross 
sections  but  no  attempt  was  made  to  include  the  auroral  component  of  absorption. 
The  measurements  of  received  clutter  power  which  were  used  were  averaged  over 
the  entire  ±  16  Hz  range  of  the  Polar  Fox  II  Doppler  processor.  Thus,  if  it  were 
desired  to  normalize  the  resultant  cross  sections  in  the  Doppler  domain,  it  would 
be  necessary  to  divide  them  by  3 2  77  ,  where  TJ  is  a  "filling  factor"  representing  the 
distribution  of  energy  within  the  range  of  the  Doppler  processor  (in  the  case  of 
uniform  Doppler  distribution,  for  example,  1}  =  1).  It  should  be  noted  also  that 
the  range  of  the  Polar  Fox  Doppler  processor  was  inadequate  to  permit  unambigu¬ 
ous  spectral  power  measurement.  Thus  the  measured  Doppler  dependence  of 
received  power  is  distorted  by  the  Doppler  fold-over  which  results  from  spectral 
power  outside  the  processor  range. 


EARTH'S  MAGNETIC  FIELD 


vs  =  r  -u-p-y 


Figure  15.  Scatter  Volume  in  the  Presence  of  Ionospheric  Refraction. 

The  effective  vertical  angle  at  the  radar  is  determined  by  the  pulse 
length 

When  the  factors  noted  in  the  preceding  paragraph  are  considered,  it  is  appar¬ 
ent  that  the  "pseudo  volumetric  cross  section"  measurements  must  be  treated  with 
caution.  Alternatively,  the  greatest  confidence  may  be  placed  in  them  in  applica¬ 
tions  which  approximately  duplicate  the  Polar  Fox  II  measurements,  especially  in 
respect  to  location  with  respect  to  the  auroral  oval  and  Doppler  processor  range. 
Since  the  present  application  fits  this  requirement,  it  is  considered  appropriate  to 
extrapolate  the  Polar  Fox  II  results,  provided  that  illuminated  volume  is  computed 
in  the  same  approximate  way,  and  auroral  absorption  is  not  included  in  the  loss 
term  of  the  radar  equation. 


The  "pseudo  volumetric  cross  sections"  from  the  KADC  Study  may  be  repre¬ 
sented  by  separate  normal  distributions  for  E-  and  l-'-layers  for  each  of  the  four 
frequencies  studied,  in  the  range  7-19  MHz  (Figures  16  to  19).  The  medians  for 
these  distributions,  at  10  MHz,  can  be  estimated  as 

~  -99  dB(m“/m\  E-region; 

-  -77.5  dB(m2/m3),F-region. 

The  frequency  dependence  of  median  o',  in  the  range  7-13  MHz,  can  be  represented 
by 

_  -3.0  dB/MHz,  E-region 
”dl  -5.  2  dB/MHz,  F-region. 


Above  about  13  MHz,  the  frequency  dependence  is  negligible. 

In  order  to  estimate  the  dependence  of  volumetric  cross  section  on  geomagnetic 

activity,  data  from  the  Prince  Albert  Radar  taken  at  400  MHz  were  used  (PARL 
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Report,  pp  7-8  and  7-9).  From  these  data  a  dependence  on  in  the  range  -  5 
can  be  inferred  by  fitting  a  straight  line  to  the  a  .  measurements; 


a  =  -113  +  2.  4  K  dBm  at  400  MHz  . 
max  p 

Using  this  same  dependence  and  assuming  that  the  Polar  Fox  II  median  measure¬ 
ments  are  representative  of  median  geomagnetic  activity,  the  above  median  "pseudo 
cross  sections"  can  be  rewritten  as 


cr'  =  -94+  2.  4  Kp  dll  (m”/m3).  E-region; 

=  -83+  2.  4  Kp  dli(m“/m\  F-region  . 

It  may  be  noted  that  Cireenwald  et  al33  found  that  the  backscattered  auroral  power 
at  College  was  approximately  proportional  to  the  square  of  the  perturbation  of  the 
horizontal  component  of  the  magnetic  field  vector,  in  qualitative  support  of  the 
above  derived  geomagnetic  activity  dependence. 


68.  Greenwald,  R.  A.,  Eckland,  \V.  I,.  ,  and  Paisley,  II.  H.  (1973)  Auroral 
currents,  irregularities,  and  luminosity,  .J.  Geophys.  Res.  78:8193. 
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T'v  {dB-nf/  m5)  ^  in-  <  (dB-rrfrm3) 


The  deviations  of  the  normal  distributions  representing  "pseudo  volumetric" 
reflectivities  may  be  taken  as 

10  dB  (E-layer), 

13  dB  (F-layer), 

more  or  less  independent  of  frequency  in  the  range  8-13  MHz. 

The  values  deduced  above  can  be  used  to  deduce  backscatter  clutter  amplitudes 
from  the  Polar  Fox  site  or  from  a  site  at  a  similar  geomagnetic  latitude.  In  order 
to  extrapolate  these  results  to  other  locations  and  times,  the  probability  of  irregu¬ 
larity  occurrence,  deduced  in  Section  4.  1  and  4.  2  can  be  used.  The  underlying 
assumption  is  that  the  volumetric  cross  section  is  proportional  to  the  critical  fre¬ 
quency  of  sporadic-E,  in  the  case  of  E-layer  clutter  and  to  the  probability  of 
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spread-F  for  F-region  clutter.  The  Prince  Albert  Radar  Study  confirms  the  first 

of  these  assumptions.  The  dependence  of  £Tmax  on  sporadic-E  blanketing  frequency 

(f_E  )  at  Fort  Churchill  is  apparent  in  that  a  increases  by  ~  22  dB  per  MHz 

increase  in  fgEg.  In  order  to  "calibrate”  the  probability  distribution  of  sporadic-E 

critical  frequency  and  relate  the  probability  that  fQEg  2  3  MHz,  used  in  Section  4.  1, 

measurements  of  f  E  were  collected  for  October  1971  from  Fort  Churchill,  Goose 
o  s 

Bay,  Ottawa,  and  St.  Johns.  The  approximate  geomagnetic  latitudes  of  these  loca¬ 
tions  are  70°,  65°,  60°  and  59°  respectively.  Figure  20  shows  the  diurnal  variation 

of  f  E  for  these  locations.  Note  that  the  median  values  for  this  month  reflect  major 
o  s 

features  of  the  space-time  distributions  of  probability  (fQEg  2  3  MHz)  in  Section  4.  1, 
particularly  the  nighttime  minimum  at  <J>  ~65°  (Goose  Bay). 

The  majority  of  E-region  clutter  was  observed  by  Polar  l-'ox  II,  at  night,  at  a 
geomagnetic  latitude  of  about  70°.  The  median  nighttime  f  E  at  Fort  Churchill 

os  gq 

<4>  -70°)  in  October  is  ~5  MHz.  Assuming  a  log-normal  distribution  for  foEg 
with  a  5  MHz  median  value,  the  3  MHz  probability  of  occurrence  would  be  -0.  83. 
This  is  very  close  to  the  value  predicted  by  the  model  described  in  Section  4.  1  for 
the  probability  that  foEg  3  3  MHz.  Thus,  a  correspondence  can  be  made  between 
the  probabilitv-of-occurrence  modelled  values  and  pseudo  volumetric  cross  section 
values.  Taking  into  account  the  statistical  distribution  of  the  Polar  Fox  11  cross- 
section  data  and  their  frequency  dependence,  the  probability  contours  of  Section  4.  1 
correspond  to  the  probability  that 

a’  2  -89dB(m2/m3)  at  10  MHz. 

V 


69.  Cherneyshev,  O.  V'.  (1968)  f  E  Distribution  curves.  Geomag,  and  A  cron. 
8:902.  °  s 
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The  correspondence  thus  obtained  between  the  probability  of  sporadic-lv  occurrence 
and  volumetric  cross  section  permits  the  calculation  of  K- region  clutter  amplitudes 
to  be  performed.  This  will  be  pursued  further  in  Sections  4.  (i  and  (i.4. 
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Figure  20.  Diurnal  Variation  of  ft)Fs  for  October  1072  for  Fort  Churchill  (Oil). 
Goose  Ha.v  (tin),  Ottawa  (Oi  l,  and  St.  Johns  <S.I> 

The  quantitative  relationship  between  the  probability  of  spread- F  (as  derived 
in  Section  4.  2)  and  median  F- region  volumetric  cross  section  is  more  diflicult  to 
establish  because  spread- 1  is  not  as  precisely  defined  as  f  Hy  considering  an 

average  ionosonde  as  a  radar  and  solving  tin*  radar  equation,  it  is  possible  to 
establish  the  minumum  detectable  F- region  clutter  signal  and  thus  the  threshold  of 
spread-F  measurement.  It  was  decided  not  to  rely  on  this  approach  because  of 
many  practical  uncertainties  and  the  spread-1  probability  contours  derived  m 
Section  4.  2  were  related  to  <j'^  median  values  bv  simple  comparison  of  the  predicted 
probability  values  and  observed  cross-section  values  for  the  1’olar  Fox  II  coverage 

area  in  October.  This  resulted  in  the  approximation  of  the  !'0  percent  probabilitv 

•  >  ) 

contours  in  Section  4.  2  bv  contours  of  median  u'  ,  which  is  -7  i.  5  dll  <m“  m  at 


10  M 1 1  /.  lor  average  K  .  Thus  tile  contour  plots  in  Section  4.  2  cun  be  used  to  com¬ 
pute  F- region  clutter  signal  amplitude  if  it  is  assumed  that  they  represent  the 
probability  that  <j'  >  -94  dll  (m“/m3)  at  10  Mllz  tsee  Appendix  D). 

As  an  example  of  the  foregoing  arguments,  consider  the  results  of  Sections  4.  1 
and  4.  2  at  <t>  =  70°  (COB)  in  June  at  00  UT,  lor  -  2.  From  the  contour  plots 

p(B)  =  0.  5  . 

p(  F)  =  0.  8  . 

Thus,  at  10  Mllz, 

prob(a^(H>  >  -H9  dB(m2/m3))  =  0.  5  . 

and 

prob(0 '(!•')  >  -94  dB(m2/m3))  =  0,8  . 

Therefore,  using  the  statistical  distributions  for  u'r  (Bland  a '(F),  their  median 
values  are 

-89  dB(m2/m3),  B- region. 


and 


-85  dB  (m2/m3),  F  -  region  , 

The  values  for  median  deduced  herein  may  be  compared  with  I'llF  results, 
in  the  case  of  B-region  clutter  at  least,  bv  using  a  crude  frequency  scaling  law 
published  by  Chestnut  et  al,  in  which  a  increases  by  ~  15  dB/400  MHz  in  the 
LI1F-IIF  range.  Both  the  Hrince  Albert  Itadu  rand  the  Perimeter  Acquisition  Radar' 

(PAID  both  measure  a  — 90  dB(m2/m3)  at  90°  aspect  angle  at  -400  MHz.  Thus, 

v  2  3  , 

at  10  Mllz,  a  -  -75  dH(m  /m  ).  The  median  "pseudo  volumetric  cross  section  ' 

v  2  3 

for  the  B-layer  is  -  -89  dB(m  /m'  1.  When  allowance  is  made  for  the  ratio  of  the 
actual  Polar  Fox  II  illuminated  volume  to  that  assumed  by  King  and  Richards*11*  and 
when  ionospheric  focussing  is  taken  into  account,  a  factor  of  -  15  dB  must  be  added 
to  the  "pseudo"  cross  section.  Thus  the  estimate  of  -89  -t  15  -  -74  dB  (m2/m3)  is 
in  good  agreement  with  the  extrapolated  UHF  results. 

70,  Chestnut,  W.(i.  ,  Hodges,  ,J.  C.  ,  and  I.eadabrand,  R.  I„  (1968)  Auroral  Back 
Scatter  Wavelength  Dependence  Studies,  KAIX  -T  R-68-286. 
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4.4  Aspect  Sensitivity  of  Auroral  Clutter 


The  subject  of  the  magnetic  aspect  sensitivity  of  auroral  scatter  has  been  ex¬ 
tensively  treated  in  the  literature  but  a  great  deal  of  confusion  still  exists  regarding 
the  quantitative  estimates  of  the  degree  of  aspect  sensitivity.  It  appears  that  much 
of  this  confusion  may  be  traced  to  the  fact  that  the  geomagnetic  field  is  distorted  by 
the  strong  ionospheric  currents  which  accompany  the  aurora  (the  very  currents 
which  generate  the  plasma  instabilities  which  in  turn  are  responsible  for  the  radio 
aurora).  Despite  the  fact  that  highly  accurate  models  are  available  to  describe  the 
earth's  magnetic  field  of  internal  origin,  it  is  difficult,  if  not  impossible,  at  pre¬ 
sent  to  correct  these  models  for  time  varying  ionospheric  current  effects.  These 
currents  may  easily  change  the  apparent  direction  of  the  magnetic  field  in  their 
vicinity  by  several  degrees,  and  since  the  apparent  radar  cross  section  of  the  aurora 
varies  by  ~  10  dB  per  degree  of  aspect  angle  variation,  this  can  be  a  major  factor  in 

cross  section  estimation. 

71  70 

Leadabrand  et  al  and  Chestnut  et  al 1  studied  the  aspect  sensitivity  of  UHF' 

auroral  clutter  near  normal  incidence  (\p  ~90°)  and  reported  a  variation  of  ~  10  dB 
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per  degree  off-orthogonality.  Unger  et  al  using  the  Prince  Albert  Radar  at  UHF 

found  a  much  smaller  value  of  ~3-6  dB/deg  in  the  range  95°  <  if/  <  105°.  The  most 

7  2 

thorough  study  of  UHF  aspect  sensitivity  appears  to  be  that  of  Mitchell  and  Brown 
who  used  the  Perimeter  Acquisition  Radar  at  Grand  Forks,  North  Dakota  in  the 
450  MHz  band.  These  authors  noted  a  great  variability  in  measured  aspect  sensitivity, 
despite  careful  attempts  to  remove  sources  of  observational  error  and  bias,  which 
they  interpreted  as  the  result  of  geomagnetic  field  distortion,  in  the  neighborhood  of 
the  auroral  scatterers,  by  the  intense  auroral  electrojet  currents.  They  reported 
values  of  aspect  sensitivity  in  the  range  0~  12  dB/deg  for  diffuse  clutter  returns  and 
as  high  as  20  dB/deg  for  discrete  returns.  It  was  noted  as  a  significant  feature  that 
discrete  clutter  returns  usually  displayed  a  higher  degree  of  aspect  sensitivity  than 
did  diffuse  returns. 

The  observations  in  the  previous  paragraph  refer  to  E- region  clutter  at  VHF  and 
UHF.  There  is  reason  to  believe  that  E- region  irregularities  may  display  different 
aspect  sensitivity  in  the  HF  band  due  to  the  much  larger  scale  size  of  plasma  irregu¬ 
larities  contributing  to  the  scatter;  it  is  by  no  means  clear  that  the  plasma  instability 
processes  operating  in  the  centimeter  and  decameter  ranges  are  identical.  Further¬ 
more  it  may  also  reasonably  be  anticipated  that  F-region  clutter  may  exhibit  differ¬ 
ent  aspect  sensitivity  to  E-region  clutter,  again  because  different  instability 

71.  Leadabrand,  R.  L. ,  Larson,  A.G.,  and  Hodges,  J.  C.  ( 1967)  Preliminary 

results  on  the  wavelength  dependence  and  aspect  sensitivity  of  radar  auroral 
echoes  between  50  and  3000  MHz,  J,  Geophys.  Res.  £2:3844. 

72.  Mitchell,  MJ.  ,  and  Brown,  J.  L.  (1976)  PAR  Auroral  Study,  Vol.  5,  Final 

Rept.  Contract  DASG-74-C-0026,  AIMlD  Systems  Lommand. 
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mechanisms  are  almost  su  'v  operating  in  the  two  altitude  regions.  Unfortunately 
there  exist  few  reliable  direct  measurements  at  UK,  owing  to  the  combination  of 
instrumental  difficulties  and  the  uncertainties  in  taking  account  of  ionospheric  re¬ 
fraction.  Bates  and  Albee1^  reported  a  value  of  -  6  dB/dog  while  MeDiarmid ‘ 

reported  a  much  lower  value  of  1.  3  -  1.5  dU/deg  at  a  low  V  UK  frequency.  As  noted 

35 

in  Section  3.  2  the  observations  in  England  of  Baggaley  can  be  interpreted  as  in¬ 
dicative  of  an  aspect  sensitivity  of  ~  10  dB/deg  for  K-region  backscatter  at  UK. 

In  view  of  the  paucity  of  experimental  evidence  and  the  apparent  great  variabil¬ 
ity  of  aspect  sensitivity  (as  evidenced  by  the  work  of  Mitchell  and  Brown) it  seems 
that  the  best  which  can  be  done  at  present  is  to  assume  a  representative  value  of 
~6  dB/deg  for  UK  aspect  sensitivity. 


1.5  Auroral  Clutter  Amplitude 

4.5.  1  HEKKAt'TED  CASK 

The  volume  of  auroral  ionospheric  irregularities  illuminated  by  the  radar  is 
determined  by  the  azimuthal  antenna  beaimvidth,  the  radar  pulse  length,  and  the 
aspect  sensitivity  of  the  scatter  from  the  irregularities  (see  Appendix  C) 


whore 


V  =  scattering  volume, 

11  =  range  to  scattering  volume, 

Ap  =  azimuthal  beaimvidth  of  antenna. 


At  =  range  of  elevation  angles  illuminating  scatter  volume. 


I*  =  group  path, 

if/  =  magnetic  aspect  angle. 


K  =  focussing  factor. 

Several  parameters  which  enter  Eq.  (1)  -  namely,  A  t  , 


ill1 

rnf  * 


and  I  must  he 


determined  by  mathematical  ray-tracing  through  a  model  ionosphere,  in  the  pre¬ 
sence  of  a  magnetic  field,  at  the  frequency  under  consideration.  In  other  words, 
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these  parameters  are  complicated  functions  of  frequency  and  the  detailed  ionospheric 
structure  between  the  radar  and  the  clutter  volume.  As  a  result  of  performing  such 
propagation  simulation  computations,  representative  values  for  these  pa  ra meters  can 
bo  established  as  follows: 

(1)  A  t  is  determined  primarily  by  the  radar  pulse  length  and  also 
by  the  range, 

that  is,  A  t  =  f<T  ,  R) 

where  T  is  the  radar  pulse  length.  Since  A  t  is  the  range  of 
elevation  angles  over  which  orthogonality  is  achieved  with 
the  magnetic  field  at  range  R,  this  parameter  is  an  approxi¬ 
mate  measure  of  the  vertical  extent  of  the  scatter  volume  (11) 

11  ~  R  ‘  At  . 


Note  that  A(  is  also  a  function  of  the  magnetic  field  geometry 

and  thus  varies  with  radar  location; 

dp 

(2)  is  an  approximate  measure  of  the  horizontal  dimension  of 

the  scatter  volume  (1.) 


A  , 


where  the  aspect  sensitivity  of  the  scatter  function  determines 
the  range  extent  over  which  seatterers  will  contribute  to  the 
scattered  power  received  at  the  radar.  The  range  of  magnetic 
aspect  angles,  A  l// ,  is  chosen  somewhat  arbitrarily  as  1  1". 
Since  a  representative  figure  for  aspect  sensitivity,  in  the  111 
range,  is  6  dll  per  degree,  this  range  effectively  defines  the 
horizontal  dimension  of  the  scatter  volume  by  its  0  dlt  points. 

Cl)  K  is  a  measure  of  the  focussing  which  results  from  the  spatial 
variation  of  ionospheric  refractive  index  and  the  curvature  of 
the  ionosphere.  It  is  essentially  the  ratio  of  energy  density 
per  unit  area  orthogonal  to  the  ray  direction  in  the  ionosphere 
to  thi'  same  energy  density  in  free  space.  Typical  values  for 
I  near  the  apogee  of  ray  trajectory  are  10-la  illt. 


The  volume -scattering  cross  section  (cr  )  is  defined  as  the  ratio 

v 

lent  point-target  cross  section  (a)  to  the  illuminated  volume. 


a 

v 


Using  the  standard  radar  equation,  a  is  expressed  by 


G) 


(4jt) 

P 


,4 


kT  H 


T 


G7T 

x 


where 

S 

N 

L 

k 

T 


PT 


is  received  signal  to  noise  ratio, 
total  system  loss, 

-■>3 

Bolzmann's  constant  (1.38  X  10  “  joule/ ’’K), 

system  noise  temperature, 

receiver  bandwidth, 

transmitted  power, 

transmit  antenna  effective  gain, 

receive  antenna  effective  area. 


The  receive  antenna  effective  aperture  can  be  expressed,  if  desired, 
its  directive  gain.  1)R,  by  the  approximate  expression 

A  2 

R  .  ,  41000  Az 

~nr  rtR  4  -  a  0  •  so  • 


where 


A  -  is  the  radar  wavelength, 

A 0  -  is  the  azimuthal  3  <IB  beamwidth, 

A 9  -  is  the  elevation  3  dB  beamwidth, 

(A0  and  A0  are  both  expressed  in  degrees).  The  use  of  l>|.  (4)  is  a 
low  frequencies  where  it  can  reasonably  be  assumed  that  the  radar  b 
limited,  so  that  antenna  efficiency  need  not  be  considered. 

Fquations  (1),  (2),  (3),  and  (4)  may  be  used  to  define  the  volume 
section  as  measured  by  a  HF  radar  as 


of  the  equiva  - 

(2) 


(3) 


in  terms  of 


(4) 
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Alternatively,  if  the  volume  scattering  cross  section  is  known  (or  can,  for  example, 
be  extrapolated  from  another  frequency  regime,  the  received  clutter  signal  can  be 
estimated  from  Eq.  (6)  by 


S 

N  ' 


2  HP 

K<Iv  A  PtGtA£  A*  ^ 

-  - £_ 


R  L  kT  BFAO 
s 


(7) 


where 

K  ■  41000  rar  1  1«- 

4.5.2  LINE-OF-SIGHT  LIMITING  CASE 

The  computation  of  received  clutter  strength  is,  in  general,  a  complicated  pro¬ 
cess,  due  principally  to  the  computation  of  effective  scatter  volume  in  the  presence 
of  ionospheric  refraction.  A  special  case  of  considerable  importance  is  found  in 
the  limit  when  the  refraction  approaches  zero,  namely,  at  night  and  at  high  fre¬ 
quencies;  this  limit  may  be  referred  to  as  the  "line -of-s ight"  limit,  bv  analogy 
with  the  UHF  radar  case.  In  this  limit  the  main  considerations  in  the  determina¬ 
tion  of  the  scatter  volume  concern  the  geometry  of  the  propagation  path  and  the  geo¬ 
magnetic  field.  Figure  21  illustrates  the  contours  of  the  geomagnetic  aspect  angle, 
t j. ,  for  this  limit,  assuming  auroral  scatterers  at  an  altitude  of  110  km. 

The  calculation  of  effective  clutter  cross  section  is  best  performed  bv  convolv¬ 
ing  the  antenna  pattern  and  the  radar  pulse  shape  (transformed  to  a  range  function) 
with  the  volume  scatter  cross  section  expressed  as  a  function  of  geomagnetic  aspect 
angle  (in  range-azimuth).  Special  attention  is  still  required  in  the  calculation  of  the 
range  of  elevation  angles  to  be  included  in  the  scatter  volume.  For  radar  antennas 
having  a  wide  elevation  beam  pattern,  the  "elevation  component"  of  the  scatter 
volume  is  no  longer  determined,  as  it  is  in  the  refracted  case,  by  magnetic  aspect 
sensitivity.  In  many  cases  it  will  be  sufficiently  accurate  to  overlay  an  equivalent 
range-azimuth  cell  (determined  by  pulse  length  ami  azimuthal  bcamwidth)  on  the 
contours  of  Figure  21,  noting  that  the  range  dimension  of  the  cell  must  be  fore¬ 
shortened  by  a  factor  0.  983  cos  t  ,  where  t  is  the  elevation  angle.  The  range  of 
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elevation  angles  is  determined  by  the  projection  of  the  pulse  length,  in  range  units, 
on  the  thin  auroral  layer.  Figure  22  shows  this  range-elevation  function.  The 
average  value  oft//  is  estimated  from  the  graphical  overlay  and  the  effective  volume 
scatter  cross  section  computed  therefrom,  using  a  figure  of  -6  dB  per  degree  off- 
orthogonalitv. 


AZIMUTH  (DEG) 


Figure  21.  Range-azimuth  Variation  of  Geomagnetic  Aspect  Angle  at  110  km 
Altitude  for  Line-of-Sight  Propagation  From  Site  at  61°  CGL,  Using  Jensen- 
Cain  Model  of  the  Geomagnetic  Field 

The  effective  scatter  volume  is 

Vg  =  R2  •  0  •  At  •  r  •  £,  (8) 

where 

A t  is  the  elevation  angle  range, 

0  is  azimuthal  beamwidth, 

T  is  pulse  length, 

c  is  velocity  of  light. 


Figure  22.  Range -elevation  Angle  Dependence  for  Clutter  at 
Altitudes  of  110  km  and  300  km  (line-of-sight  propagation) 


Note  that  the  radial  component  of  V  may  be  limited,  in  the  case  of  very  long 
pulse  length  and/or  steep  incidence  on  the  ionosphere,  not  by  pulse  length  but  by 
the  thickness  of  the  scattering  laver  in  the  E  region,  which  is  approximately  10  km. 
Thus  if  r  •  exceeds  10  sec  8  where 

8  -  sin  1  (0.983  cos  <  )  (9) 

and  f  =  elevation  angle,  then  the  smaller  of  these  two  values  is  used  in  calculating 
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where 


P 

f 

1' 

o 

Z 


max 

Z 

Y 


m 


r 

T  ■ 

o 

operating  frequency, 
maximum  plasma  frequency, 

height  of  layer  maximum, 

altitude  measured  from  Z 

max 

layer  semi -thickness. 


After  refraction  in  this  layer,  the  ray  reaches  orthogonality  with  the  magnetic 
field,  according  to  Snell's  Law,  where 


n(Z) 


sin  j 

o 

sin  I 


(11) 


where 

:>o  =  angle  of  incidence  on  layer, 

I  magnetic  inclination. 


The  minimum  slant  range  at  which  orthogonality  is  achieved  corresponds  to 
rays  which  reach  the  peak  of  the  layer  and  are  almost  penetrating,  that  is,  Z  -  0 
in  Kq.  (10). 

Assuming  a  given  range,  U,  the  angle  of  incidence  on  t he  curved  laver  is  given 
by 


cos  0 

o 


H2  +  2RehB  ♦  hB2 

— 2TTTTT 


-  _r_fvr 

e  B 


(12) 


where 


R  =  earth  radius  (6380  km), 
e 

hB  -  height  of  underside  of  layer. 


that  is 


h„  (Z  -  Y  ). 
B  max  m 
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Solving  Kqs.  (10),  (11),  and  (12)  yields  the  slant  range  to  the  layer  bottom  for 
which  orthogonality  is  attained;  it  is  then  necessary  to  compute  the  component  of 
virtual  slant  range  for  the  refracted  part  of  the  raypath,  within  the  layer.  Let 
this  component  be  p'(Z). 

Then 

Z 

p'(Z)  =  J  ^ 

h  q 
B 


where 


2  2  .2 
q  *  n  y  sin  <[> 


that  is. 
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Upon  integrating  Fq.  (13), 
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The  virtual  slant  range  to  the  orthogonality  point  is  then  given  bv  It  •  p ' ,  from 
Eqs.  (12)  and  (14).  The  minimum  slant  range  is  the  value  of  It  *  p '  for  Z  =  0, 
that  is,  corresponding  to  rays  reaching  the  peak  of  the  layer.  This  is  plotted  as 
a  function  of  p  in  Figure  22a  for  an  assumed  layer  bottom  height  of  300  km,  semi- 
thickness  of  100  km,  and  a  magnetic  inclination  of  70’. 
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f  •  fi  'i  iaai 


km'  100  km 

I  >70* 


Figure  22a.  Minimum  V  irtual  Slant  Mange  for 
F -region  Clutter  as  a  Function  of  Frequency, 
Normalized  to  Critical  Frequency  of  a  Sphericallv 
Stratified  Caver.  Vm  -  parabolic  lavor  semi- 
thickness,  h..  height  of  laver  bottom, 

I  =  geomagnetic  inclination 


5.  DOPPI.KK  SHIFT  MODF.I. 

As  noted  in  Section  2.2,  the  Doppler  shill  model  is  based  on  the  theoretical 
29 

work  of  Mnoda,  referred  to  therein,  with  some  relatively  minor  modifications 
which  were  found  necessary  in  order  to  being  Mnoda' s  velocity  predictions  into 
better  agreement  with  available  observations.  iWaeda's  calculated  values  a  re 
assumed  to  be  appropriate  for-  a  planetary  magnetic  index  — K  -  2,  that  is,  an  aver 
age  value.  Values  of  the  northward  and  eastward  components  of  plasma  \cloeits 
(V  ^  and  V  j  )  in  tile  CCS  were  tabulated  at  intervals  of  geomagnetic  latitude  ,ud  d 
2-hr  intervals  of  geomagnetic  time.  A  numerical  interpolation  procedure  was  then 
used  to  compute  these  velocity  components  at  arbitrary  latitude  owl  turn  .  V  sepa¬ 
rate  model  was  employed  for  li- region  and  F- region  velocities  respective  l\.  II  w  a 
assumed  that  the  entire  velocity  pattern  expanded  lat  itudiilallv  b\  2  per  amt  of  , 
reflecting  the  expansion  of  the  auroral  oval  with  increasing  magnet  ii  activitv. 
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Table  1.  E-Kegion  l'lasnia  Velocities  in  CCIS 


a 

O 

O 

CO 

II 

•e- 

<t>=  75° 

<t>=  7 

0° 

<t>=  65° 

$  =  60° 

VN 

Q 

VN 

VE 

VN 

B9 

VN 

E9 

VN 

ra 

0 

-117 

-305 

-229 

-179 

-171 

-141 

-  70 

-  83 

0 

50 

2 

115 

-226 

-149 

-229 

-260 

-  80 

-211 

142 

-38 

32 

4 

197 

-160 

-  22 

-254 

-340 

173 

-192 

167 

-41 

29 

6 

244 

-  70 

106 

72 

-352 

365 

-135 

122 

-41 

29 

8 

182 

0 

57 

248 

-199 

264 

-  75 

51 

-41 

29 

10 

173 

56 

-  94 

352 

-  63 

173 

-  66 

31 

-50 

0 

12 

224 

156 

27 

216 

-  43 

59 

-  34 

-  43 

-32 

-38 

14 

194 

240 

182 

0 

42 

-  26 

40 

-  61 

43 

-25 

16 

-  60 

341 

220 

-  8b 

287 

-  166 

58 

-  70 

32 

-38 

18 

-135 

318 

218 

0 

573 

-209 

98 

-108 

38 

-32 

20 

-351 

88 

-128 

0 

427 

-288 

240 

-  87 

41 

29 

22 

-250 

-210 

-179 

-  32 

95 

-224 

205 

-119 

49 

9 

5.2  E-Region  V»locity  Model 

The  model  for  F-region  plasma  velocity  is  similar  to  that  for  the  E- region, 
except  that  the  latitudinal  asymmetry  associated  with  the  auroral  oval  is  taken  into 
account  explicitly. 

Let  a  geomagnetic  latitude  boundary  be  defined  as  a  function  of  geomagnetic 
time  by 


=  72  -  K  -  5 
F  p 

(1)  When  $  =  4>j,. 


cos 


7! 

12 


T  . 


VN  =  400  cos  jj-  (T  -  18);  6  <  T  <  18 

=  400  cos  |  (T  -  18);  T  2  ia;  T  6  . 

VE  =  0 

(2)  When  <I>  =  -  3 

vN  -  o 

V.,  =  500  +  500  cos  £  (T  -  8);  0  <  T  <  12 

li.  b 

=  -500  -  500  cos  |  (T  -  HI);  12  '  T  ’  2-1  . 
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(3)  When  4  =  4-F  -  7 


N  ' 

0; 

T  4  12  and  4  00 

E  = 

200  +  70  cos  |  (T  -  6); 

0  <  T  <  12 

= 

-200  -  70  cos  |  (T  -  18); 

1 2  <  T  <  24 

For  T  =  12,  VM  =  130,  V„  =  0 
N  E 

For  T  =  00,  VN  =  -130,  VE  =  0  . 

(4)  When  4=  ®F  -  12 

VN  =  0:  T  4  12  and  4  00 

VE  =  60  cos  |  (T  -  6);  0  <  T  <  12 

=  -60  cos  |  (T  -  18)  12  <  T  <  24 

For  T  =  12.  V’  =  60,  V„  =  0 

iN  L 

For  T  =  00.  VN  =  -60.  Vj,  =  0  , 

(5)  For  4  >85° 

VN  =  -400  cos 

VE  =  400  cos  j2  (T  -  18)  . 

(6)  For  4=  80°,  V  and  \’E  are  given  by  the  values  in  Table  2. 

As  was  the  ease  for  E- region  velocity,  a  computer  program  was  written  to 
estimate  the  velocity  for  F-region  plasma  at  arbitrary  latitude  and  time,  using  the 
analytic  relationships  noted  above  and  performing  numerical  interpolation  where 
the  velocity  components  are  specified  numerically.  The  only  magnetic  activity 
dependence  contained  in  both  the  E-  and  I- -region  models  is  the  latitudinal  expan¬ 
sion  effect.  Although  various  reports  have  appeared  in  the  literature  regarding  a 
dependence  of  velocity  on  magnetic  activity,  this  must  still  be  regarded  as  an 
unresolved  issue.  The  careful  study  by  l  nger  et  al'  ,  using  the  backscatter  radar 
technique,  did  not  reveal  any  pronounced  systematic  dependence  of  Doppler  shift  on 
geomagnetic  activity.  Although  it  is  clear  that  large  and  often  erratic  variations 
in  plasma  velocity  do  occur,  the  model  described  here  is  regarded  as  a  reasonable 
characterization  of  average  velocities. 
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Tabli*  2.  1'- Region  Plasma 

Velocities  at  <t>  =  110° 


T 

A 

S: 

0 

-400 

0 

-> 

-  2<i  7 

-297 

4 

-  35 

-398 

0 

-400 

8 

104 

-386 

10 

160 

-362 

12 

400 

0 

14 

16!) 

36  2 

16 

104 

386 

la 

0 

400 

20 

-  35 

398 

22 

-26  7 

2J»V 

a. 3  linpftlt-r  Sliifl  (li>i»|»iliil ions 

The  velocity  model  described  in  Sections  -1.  I  and  ■!.  2  gnes  tin-  north  ,nd  cast 
components  of  K-  and  1  -region  plasma  motion  in  the  t  ils,  m  order  to  compute  the 
resultant  Doppler  effect  on  the  radar  signal,  it  is  iieeessarv  to  perform  appropriate 
coordinate  transformations  So  as  first  to  calculate  the  radial  component  of  motion 
relative  to  tile  radar,  figure  21  illustrates  the  situation  tin  the  approximation  of 
plane  geometry).  M  and  (1  are  n-spee-1 1\ els  the  geomagnetic  md  geographic  poles 
and  1‘  is  the  location  of  the  clutter-producing  plasm,,,  s  gs  the  loe  ,tion  of  the  radar 
site.  \  and  \  | .  a  re  the  velocitv  component  -  p  im  id,  d  t>\  1  he  mode  I,  whose  result¬ 
ant,  \  ,  is  the  actual  plasma  \eloe  lU  vector.  1  he  components  of  \  ill  the  geodetic 
system  are  \  ^  and  \  j  and  the  radial  veloeit\  ,  oinpotic  nt  rel  dive-  to  the  radar  is 
\  ^.  It  is  roqui  red  to  dote  run  lie  \  ^  in  te  rms  ot  \  \  ,  ,nd  hot  h  the  g  cog  rapine 

and  geomagnetic  locations  of  both  1’  md  s. 

As  noted  m  the  preceding  paragraph,  the  plain-  geometric  t  r<  atmnit  of  this 
problem  is  iilaecur-de  and  is  mentioned  oiiD  tor  purposes  ot  illustration.  In  fact, 

.1  spile  lie  al  f  ingonoiliet  lie  solution  Is  lie,  ISVU'C  Cue  \ppendlN  li).  \  computer 
prog  r- 1 m  was  w  rittcii  to  pe  r!o rn ,  t  In  neercs  ,  r  oionlin  ,le  1 1  uistnrn ,  at  ions  and 
coupled  to  the-  program  generating  tin-  \  clot  it .  component-.  I  lie  re  sultant  radial 
component  of  \  e  -  loe  It  V  rel  at  U  e  t  ”  t  lie  a  It  e  ,  1  lit.  i  \  ;  le  ,  .  i  A  1  W  .1  a  e' cm  put  ed  .  e 
function  ol  range-  end  ,/imuth  aul  gieuu  egneti*  otr.itc,  torteotn  |,-  ,m|  I  -regions. 
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Figure  24.  Computation  of 
Kadial  Component  of  l’lasma  Velocity 
for  Site  at  S,  When  Velocity 
Components  (Vj^,  V'^)  are 
Specified  in  Geomagnetic  Coordinates. 
M  is  geomagnetic  pole,  G  is 
geographic  pole.  Plane  geometry 
is  assumed 


The  radial  component  of  velocity  was  transformed  to  a  Doppler  shift,  normalized 
to  the  radar  frequency  (units  of  Hz /MHz),  and  using  the  formula 


where 


f  =  radar  frequency, 

Vg  =  radial  velocity  component, 

C  =  velocity  of  light. 

Note  that,  by  convention,  a  negative  Doppler  shift  signifies  an  approaching  velocity. 

Figures  25  and  26  show  typical  results  of  these  computations  for  F-  and  F- 
region  irregularities  respectively  (for  a  local  meridian  time  at  the  site  of  05  1U1.T). 

It  is  apparent  that  much  complex  structure  exists  in  these  Doppler  maps,  and  the 
rough  outline  of  the  auroral  oval  can  be  seen  clearly  in  the  !•  -region  map.  Also 
noteworthy  is  the  range  of  Doppler  frequencies  predicted  by  the  model.  It  is  by  no 
means  unusual  to  see  Doppler  offsets  of  5  Hz/MHz,  and  values  as  high  as 
6-7  Hz/MHz  are  occasionally  present  in  localized  "hot-spots."  Although  computa¬ 
tions  were  performed  for  two  levels  of  geomagnetic  activity  (K^  =  2,4)  only  the  re¬ 
sults  for  average  activity  (K^  =  2)  are  presented  here.  Hecause  of  the  manner  in 
which  geomagnetic  activity  enters  the  model,  the  effect  of  increasing  is  primarily 
to  expand  the  auroral  oval  by  1.  5°  per  unit  K  .  A  qualitative  estimate  of  the  effect 
of  this  on  the  Doppler  maps  can  be  made  bv  noting  the  position  of  the  auroral  oval 
and  its  related  Doppler  features  and  applying  an  appropriate  "expansion  factor". 


Translation  of  these  maps  to  actual  radar  frequencies  is  a  simple  linear  scaling 
procedure,  based  on  the  defining  equation. 


Figure  25.  Contours  of  E-region  Doppler  Shift  Computed  for  Site  at  (44.  5“  N, 
67.  5°  W)  at  0730  UT.  Contour  units  are  Hz/MHz  (that  is.  normalized  to  radar 
frequency).  Azimuth  is  shown  in  10°  intervals  with  the  horizontal  line  indi¬ 
cating  geographic  East -West.  Kange  is  shown  in  5°  (that  is,  556  km)  intervals 


Figure  26.  Same  as  Figure  25,  Except  for  F-region  Clutter 


As  an  example  of  the  use  of  the  Doppler  model.  Figure  27  shows  the  variation 

with  Universal  Time  of  Doppler  shift  at  10  MHz,  for  E -region  clutter  at  a  range  of 

600  km  from  a  site  located  in  Northern  Scotland.  The  dashed  curve  shows  the  time 

dependence  for  a  northward -looking  radar,  while  the  solid  curve  shows  the  similar 

dependence  for  a  look-azimuth  of  -60°.  Note  the  characteristic  reversals  of  Doppler 

direction  around  noon  and  midnight  respectively.  This  has  been  observed  in  many 
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experimental  studies.  The  results  of  Leadabrand  et  al  at  Fraserburgh,  Scotland 
and  of  Unwin'"  are  in  reasonable  quantitative  agreement  with  the  predictions  de¬ 
veloped  here.  Figure  28  shows  similar  calculations,  based  on  the  model,  for  the 
F-region  clutter.  Note  the  generally  larger  values  of  Doppler  shift,  the  greater 
azimuthal  difference,  and  again  the  diurnal  pattern  exhibiting  reversals  near  noon 
and  midnight. 

A  further  example  of  the  use  of  the  Doppler  models  is  illustrated  in  Figures 
29  and  20,  for  E-region  and  F-region  clutter  respectively.  These  figures  show  the 
azimuthal  variation  of  Doppler  shift  at  00  UT  and  03  UT  at  10  IvlHz  for  a  range  of 
600  km.  Note  that  E-region  velocity  at  midnight  reverses  close  to  the  direction  of 
magnetic  north  as  found  in  several  experimental  studies. 

5.4  Doppler  Spread 

The  spectrum  of  the  signal  backscattered  by  auroral  irregularities  indicates 
a  kind  of  turbulent  motion  of  the  scatterers  with  a  random  spread  of  Doppler  fre¬ 
quencies  which  can  often  be  characterized  bv  a  simple  statistical  distribution  func¬ 
tion.  Statistical  data  concerning  the  spectral  distribution  of  HF  signals  back- 
scattered  by  the  aurora  are  relatively  rare.  This  is  accentuated  bv  the  fact  that  the 
present  application  requires  that  the  spectrum  be  measured  over  a  very  large 
dynamic  range  (up  to  60  dB  or  more).  Very  few  experiments  have  been  conducted 
with  equipment  capable  of  achieving  such  large  dynamic  range.  The  Polar  Fox  II 
experiment,  which  did  possess  a  large  dynamic  range,  included  measurement  of 
the  spectral  width  at  a  number  of  amplitude  levels  below  the  peak  return  near  zero 
Doppler  shift.  Unfortunatelv,  this  data  has  never  been  analyzed  in  a  form  suitable 
for  the  present  purpose.  A  study  of  Doppler  spread  data,  using  original  Polar 
Fox  II  data  tapes  has  been  initiated  at  RADC,  but  no  results  are  available  for  inclu¬ 
sion  in  this  report. 

75.  Leadabrand,  R.  L.  ,  Schlobohm,  J.C.  ,  and  Baron,  M..J.  ( 1965)  Simultaneous 
very  high  frequency  and  ultra  high  frequency  observations  of  the  aurora  at 
Fraserburgh,  Scotland.  .1.  Geophys.  Res.  70:4235. 
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Figure  29.  Azimuthal  Dependence  of  F -region  Doppler  Shift  at  10  MHz,  for 
Universal  Times  00  and  03,  From  Site  at  61"  CGL  Constant  600  km  Range 


Figure  30,  Same  as  Figure  29,  Kxeept  for  F -region 
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The  Polar  Cap  III  experiment  also  possessed  sufficient  dynamic  range  to  per¬ 
mit  measurement  of  spectral  spread  down  to  levels  of  the  order  of  -60  dB  below 
ground  clutter.  Again,  however,  the  documented  results  of  this  experiment  do  not 
permit  a  useful  model  to  be  constructed  which  meets  the  present  requirements. 

Both  Polar  Fox  II  and  Polar  Cap  III  measured  the  spectral  properties  of  HF  signals 
forward  scattered  over  the  skywave  path.  It  is  considered  risky  to  apply  these 
results  to  auroral  backscatter,  due  to  the  important  differences  in  geometry  and 
to  the  greatly  different  signal  levels  involved. 

In  the  absence  of  relevant  HF  data  on  spectral  spread  of  auroral  backscatter, 

it  has  been  decided  to  adopt  a  preliminary  model  based  on  UHF  radar  data,  speci- 
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fically  data  from  the  Chatanika  radar  operated  by  SRI.  This  model  can  be  re¬ 
placed  by  another,  based  on  Polar  Fox  II  data,  when  the  appropriate  analysis  has 
been  completed  at  RADC.  Using  spectral  spread  data  from  the  above  reference,  the 
width  of  the  Doppler  spectrum  of  diffuse  auroral  backscatter  may  be  conveniently 
approximated  by  a  Gaussian  function  having  the  form 

AV  =  750  exp  -  (V2/a2) 

where 

AV  =  velocity  spread  at  half-power  points  of  the  velocity 
spectrum  (m/sec), 

V  =  velocity  at  which  the  backscattered  power  is  maximum  (m/sec), 
a  =  300  m/sec. 

This  functional  form  reflects  the  observations  that  the  width  of  the  Doppler  spectrum 
decreases  as  the  spectral  peak  moves  away  from  zero  Doppler  shift. 

6.  PREDICTABILITY  OF  AURORAL  CLUTTER 

Since  auroral  radar  clutter  is  intimately  related  to  electric  fields  and  currents 

in  the  magnetosphere  — ionosphere  system,  the  question  of  its  predictability  hinges 

largely  on  the  predictability  of  these  phenomena.  It  is  well  known  that  the  electric 

fields  and  convection  of  plasma  in  the  magnetosphere  are  in  turn  dependent  upon 

conditions  in  the  interplanetary  medium  and  ultimately  on  the  degree  and  nature  of 

activity  on  the  sun,  so  that  the  quet  'on  of  predictability  is  ultimately  concerned 

with  the  extent  to  which  solar  disturbances  can  be  predicted.  Although  consider  - 

able  progress  has  been  made  in  understanding  solar-terrestrial  relationships  in 

7  6.  Moorcroft,  D.  R.  ,  and  Tsunoda,  H.  T.  ( 1978)  Rapid  scan  doppler  velocity  maps 
of  the  UHF  diffuse  and  radar  aurora,  J.  Geophys.  Res.  83:1482. 

- > -  ^ ^ 


recent  years,  the  extent  of  knowledge  on  this  subject  falls  far  short  of  a  capability 
for  long  term  prediction  of  auroral  clutter  activity,  with  onlv  a  few  well-defined 
exceptions.  At  present,  and  in  the  forseeable  future,  the  best  that  c  an  be  hoped 
for  is  prediction  on  a  scale  of  ~  1  hr,  based  on  appropriate  monitoring  of  condi  - 
tions  in  the  near-earth  environment.  A  notable  exception  to  this  conclusion  is  the 
fact  that  large  scale  geomagnetic  storms  exhibit  certain  repeatable  properties  which 
can  be  used  for  prediction  purposes. 

Global  geomagnetic  storms  are  associated  with  a  characteristic  class  of  dis¬ 
turbance  in  the  solar  corona  and  chromosphere.  They  often,  but  not  always  follow 
another  kind  of  solar -geophysical  event  known  as  a  Polar  Cap  Absorption  (PCA) 
event.  PCA's  are  the  result  of  energetic  particles  (mostly  protons)  of  solar  origin 
which  enter  the  earth's  atmosphere  in  the  polar  cap  regions  and  produce  enhanced 
D-region  ionization— hence  their  name.  They  are  usually  followed,  after  an  interval 
of  one  to  several  days,  by  a  world-wide  geomagnetic  storm  associated  with  (among 
other  things)  enhanced  auroral  clutter.  These  storms  may  last  from  one  to  ten  davs, 
the  latter  being  an  extreme  case.  An  interesting  feature  of  this  sequence  of  events 
is  that  the  solar  active  regions  which  are  initially  responsible  for  the  terrestrial 
events  are  often  long-lived.  If  they  persist  for  a  time  longer  than  the  solar  rota¬ 
tional  period  (27  days)  they  have  a  high  probability  of  generating  another  geomag¬ 
netic  storm  on  successive  solar  rotations.  In  fact  it  is  quite  common  for  these 
solar  active  regions  to  persist  for  several  solar  rotations,  giving  rise  to  recurrent 
geomagnetic  storms  at  27 -day  intervals,  although  generally  of  decreasing  intensity 
as  the  solar  active  region  weakens  and  disperses. 

6.1  Statistical  Considerations 

Apart  from  the  recurrent  world-wide  geomagnetic  storms,  which  are  in  tiny 
case  relatively  infrequent  events,  the  major  problem  of  predictability  of  auroral 
clutter  is  associated  with  the  predictability  of  auroral  substorms  — a  much  more  fre¬ 
quent  event,  confined  to  the  vicinity  of  the  auroral  oval.  The  most  readily  available 
indicator  of  substorms  is  the  disturbance  of  the  geomagnetic  field  recorded  bv  the 
ground  based  magnetometers.  Various  indices  of  geomagnetic  disturbances  have 
been  developed;  the  Auroral  Electrojet  (AE)  index  is  derived  from  selected  high 
latitude  observatories  while  the  planetary  geomagnetic  index  (Kp)  is  intended  to  be 
representative  of  a  more  global  level  of  magnetic  activity.  The  statistics  of  Kp 
have  been  studied,  in  recent  years,  by  Elkins*®  and  Cage  and  Zawalick.  77  Pazich7® 


77.  Cage,  A,  L.  ,  and  Zawalick,  E.  J.  (1972)  A  Discussion  of  the  Geomagnetic  In- 

dices  Kp  and  A  .  1932  to  1971,  AKCRt  .-1 2- 0693,  AD  756828. 

78.  Pazich,  P7M.  ( 1&I6)  Conditional  Probabilities  of  the  Geomagnetic  Index  A  , 

AFGL-TR-76-0034,  AD  A023675.  B 
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has  studied  the  statistics  of  A  (an  average  daily  geomagnetic  index  related  to  K  ). 

79  P  P 

Elkins  has  studied  statistics  of  AE.  The  results  of  all  these  studies  can  be 

applied  to  the  statistical  determination  of  the  likelihood  of  auroral  geomagnetic 

disturbance  and  thereby,  through  the  use  of  the  model  developed  in  Section  4  for 

example,  to  the  statistical  prediction  of  auroral  clutter. 
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Elkins1  showed  that  the  probability  distribution  of  over  a  38-year  period 
was  approximately  Gaussian,  with  a  median  value  of  1.  5,  and  upper  and  lower 
deciles  of  3.7  and  0  respectively  (in  the  range  -  6).  There  are,  however,  pro¬ 
nounced  solar  cycle  and  seasonal  dependences  in  magnetic  activity  as  shown  by 
77  78  79 

Cage  and  Zawalick,  by  Pazich,  and  by  Elkins  .  The  maximum  in  the  solar 
cycle  trend  of  geomagnetic  activity  follows  the  peak  sunspot  number  by  about  1-2 
years  while  the  seasonal  effect  is  characterized  by  maxima  in  the  equinoxes.  The 
variation  of  average  yearly  over  a  typical  solar  cycle  ranges  from  ~  1.  6  at 
~  1-2  years  following  solar  minimum  to  ~2.8  following  solar  maximum.  The  aver¬ 
age  seasonal  variation  of  over  a  40-year  period  shows  maxima  of  about  2.  5  in 

March  and  September  with  minima  at  '•2.0  in  winter  and  2.  1  in  summer. 

16  79 

Elkins  '  3  has  used  the  statistics  of  to  estimate  the  number  of  periods  in 
a  10-year  interval  for  which  remains  above  a  particular  value.  His  results  are 
summarized  in  Table  3,  where  the  number  of  periods  of  duration  (T)  hours  for  which 
Kp  5  (Kp)o  is  tabulated.  This  data  can  be  used  to  establish  the  frequency  with  which 
a  radar  might  be  expected  to  experience  severe  clutter,  in  terms  appropriate  to 
operational  analysis. 


Table  3.  Duration  of  Magnetically  Disturbed  Periods 


T(hrs)/(Kp)0  2 

3 

4 

5 

6 

3 

11.  000 

5600 

2200 

800 

300 

6 

1800 

620 

210 

63 

18 

9 

300 

90 

26 

5.  2 

1.  3 

12 

70 

16 

4.  5 

0.  7 

0.  1 

6.2  Interplanetary  Parameters 

The  configuration  of  the  magnetosphere  is  known  to  depend  strongly  on  the 
direction  of  the  magnetic  field  in  the  interplanetary  medium.  It  has  been  found, 
for  example,  that  the  equatorial  boundary  of  the  auroral  oval  moves  equatorward 
by  an  average  of  0.  5-0.  6  degrees  per  gamma  increase  in  the  southward  component 


79.  Elkins,  T.J.  (1973)  An  Empirical  Model  of  the  Polar  Ionosphere,  AKCRE-TR- 
73-0331,  AD  76624TT 
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of  the  interplanetary  magnetic  field— IMF.  '  Many  authors  have  noted  the 

association  between  the  southward  turning  of  the  IMF  and  the  occurrence  of  auroral 

substorms.  Arnoldy*^2  showed,  for  example,  that  a  southward -turning  IMF  often 

precedes  substorms  by  about  one  hour.  Estimates  of  the  interval  between  an  IMF 

event  and  its  related  auroral  oval  manifestations  vary,  but  are  usually  less  than 
83 

one  hour.  Pike  et  al  for  example,  reported  that  the  Explorer  35  Satellite,  at 
geocentric  distances  of  35-60  earth  radii,  detected  southward  IMF  turning  followed 
be  equatorward  shift  of  the  equatorward  oval  boundary  after  a  time  interval  of 
10-30  minutes.  It  therefore  seems  that  with  an  interplanetary  satellite  at  reason¬ 
able  distance,  provided  that  it  is  suitably  located  relative  to  the  sun,  a  reasonably 
confident  prediction  of  auroral  activity  could  be  made  with  10-60  min  typical  ad¬ 
vance  warning. 
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Recently,  Akasofu  has  claimed  to  have  found  a  quantitative  relationship  be¬ 
tween  the  auroral  eiectrojet  magnetic  index  and  parameters  of  the  interplanetary 
medium,  as  expressed  by 

AE  -  BV2  sin4  61 2 

where 

AE  =  auroral  electrojet  index, 

B  =  IMF  intensity, 

V  =  solar  wind  velocity, 

6  =  direction  of  IMF  relative  to  the  ecliptic  plane. 

If  confirmed  by  subsequent  investigation,  this  relationship  would  be  useful  in  pre¬ 
dicting  not  only  the  occurrence  but  also  the  intensity  of  auroral  clutter. 

80.  Kamide,  Y.  ,  Burch,  J.  L.  ,  Winningham,  .1.  D.  ,  and  Akasofu,  S-l  (1976) 

Dependence  of  the  latitude  of  the  cleft  on  the  interplanetary  magnetic  field 
and  substorm  activity,  ,).  Geophys.  Res.  81:698. 
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81.  Kimide,  Y.  ,  and  Winningham,  .1.  D.  (1977)  A  statistical  studv  of  the  instantan¬ 

eous  nightside  aurural  oval:  The  equatorward  boundary  of  electron  precipita¬ 
tion  as  observed  by  the  Isis  1  and  2  satellites,  .1.  Geophys.  Res.  82:5573. 
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82.  Arnoldy,  R.  L.  (197  1)  Signature  in  the  interplanetary  medium  for  substorms, 

J.  Geophys.  Res.  76:S189, 
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6.3  Satellite  Imagery 

The  advent  of  large  scale  auroral  imaging  from  a  satellite  platform  has  not 
only  increased  understanding  of  the  dynamics  of  the  aurora,  but  also  provides  a 
basis  for  short-term  auroral  activity  prediction.  The  U.  S.  Air  Force  operates 
the  Defense  Meteorological  Satellite  Program  (DMSP)  in  which  two  satellites,  in 
approximately  circular  1000  km  orbits,  record  images  of  large  portions  of  the 
north  and  south  polar  regions  on  a  routine  basis.  These  satellites  have  an  orbital 
period  of  '90  min  and  are  9CF  apart  in  longitude,  that  is,  6  hr  apart  in  local  time, 
one  being  in  the  noon -midnight  meridian  and  the  other  in  the  dawn-dusk  time  zone. 
Figure  31  shows  a  typical  DMSP  image,  recorded  over  northern  Kurope  at 
'0120  UT  on  12  February  1975.  The  outlines  of  the  U.  K.  and  of  the  Furopean 
continent  are  clearly  visible  in  the  city  lights.  A  moderately  active  aurora  is 
visible  to  the  north  of  the  U.  K.  The  use  of  this  kind  of  real-time  data  lies  in  the 
known  fact  that  auroral  substorms  lend  to  commence  in  the  midnight  sector  of  the 
auroral  oval  and  subsequently  propagate  around  the  oval  in  a  time  interval  '  1  hour. 
Thus  a  satellite  ineasurement  in  the  midnight  sector  of  the  auroral  oval,  at  the 

time  of  substorm  onset,  can  provide  some  warning  of  subsequent  auroral  activity 
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at  other  longitudes  within  the  next  30  min  to  1  hour.  Akasofu  '  has  published  a 
detailed  description  of  the  development  of  auroral  and  polar  magnetic  substorms. 
While  much  apparently  remains  to  be  learned  about  these  complex  phenomena,  the 
present  level  of  understanding  is  sufficient  to  permit  a  short  term  predictive  capa¬ 
bility  to  be  developed,  based  on  suitable  satellite  observations  in  real  time. 

6.4  (.round- Reid  Magm-lomrlers 

Suitably  located  magnetometers  can  be  used  to  indicate  the  location  and  inten¬ 
sity  of  the  auroral  electrojet  and  thence  the  location  and  intensity  of  auroral  clutter. 
The  use  of  such  magnetometer  data  in  a  predictive  mode  is  subject  to  similar 
limitations  as  those  relating  to  satellite  observations.  In  fact,  since  a  ground- 
based  magnetometer  measurement  is  a  relatively  localized  one  compared,  for 
example,  to  DMSP  imagery,  its  usefulness  mav  be  quite  limited  unless  a 
meridional  chain  of  magnetometers  was  employed.  Meridional  chains  of 
magnetometers  have  been  established  in  Alaska  and  North  America  for  research 
purposes. 

85.  Akasofu,  S-I.  (  1968)  Polar  and  (ieomagnetic  Substorms,  Springer -\ erlag, 

New  York.  . 


Figure  31.  DMSP  Photograph  o!'  Northern  Kurope 
Showing  lUoderatelv  Active  Aurora 

6„'i  backseat  ter  Sounder  in  I  Ilf  I'nlar  Cap 

One  of  the  most  promising  techniques  tor  continuous  monitoring  ol  auroral  _ 
clutter  is  the  location  of  a  HI-  backscatter  sounder  close  to  the  geomagnetic  pole. 
Because  of  the  long  range  at  which  both  K-  and  K-region  backscatter  can  be  ob¬ 
served  at  HP.  due  to  ionospheric  refraction,  the  entire  auroral  oval  can  be  moni- 

i*  C 

tored  from  a  single  location  (Bates  ct  all.  Similar  remarks  applv  to  tin-  question 
of  predictability,  based  on  this  kind  of  data,  as  for  other  techniques.  1  lie  prob¬ 
lems  of  data  reduction  and  interpretation  with  this  technique  are  likelv  to  I..-  more 
severe  than  with  satellite  or  magnetometer  data. 

86.  Bates,  M.  p.  ,  Akasofu,  S-l.  .  and  Kimball,  D.S.  (P'T.ll  first  results  I  mm 
the  north  polar  auroral  radar,  .1.  Cieophvs,  lies.  7H;.t85r. 
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Appendix  A 

F-Region  Irregularities 

9  1  OO 

The  "permanent  maximum"  in  spread-F  noted  by  Penndorf  and  Tao  appears 

to  be  an  artifact  caused  by  observational  bias.  The  region  of  intense  F'-region 

irregularities  in  the  dayside  auroral  oval  is  known  to  be  caused  by  direct  injection 

into  the  earth's  atmosphere  of  thermalized  interplanetary  plasma,  with  peak  flux 

6  3 

in  the  10 ~  100  eV  region.  Whalen  and  Pike  have  shown  that  this  region  of  irregu¬ 
larities  extends  for  at  least  several  hours  on  either  side  of  magnetic  midday.  Direct 
satellite  measurements  of  the  particle  fluxes  which  produce  the  irregularities  have 
been  reported  by  Maehlum87  and  Burch.  64  Both  authors  show  that  there  are  marked 
variations  in  the  flux  with  the  tilt  angle  between  the  geomagnetic  axis  and  the 
ecliptic  plane.  Specifically,  Maehlum  shows  a  pronounced  variation  with  Universal 
Time  of  the  high  latitude  daytime  boundary  of  precipitation  between  -74°  and  -84° 
in  summer,  while  Burch  shows  a  smaller  but  significant  variation  (about  4°  over 
the  course  of  a  year)  for  the  low  latitude  boundary.  The  highest  dayside  location  of 
the  poleward  boundary,  observed  by  Maehlum,  occurs  at  18  UT,  while  the  lowest 
occurs  at  8  UT,  Thus  there  is  a  variation  with  Universal  Time  of  the  latitudinal 
extent  of  the  precipitation  region  of  -8°  -  10°,  with  the  maximum  width  at  18  UT. 
Maximum  spread-F  should  thus  be  observed  at  the  location  where  18  UT  corres¬ 
ponds  to  magnetic  local  noon,  namely  the  Foxe  Basin  in  Northern  Canada.  At  this 
location  the  irregularity  region  is  most  extensive. 


87,  Maehlum,  B.  N.  (1968)  Universal-time  control  of  the  low-energy  electron 
fluxes  in  the  polar  regions,  J,  Geophys.  lies.  73:3459, 
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Appendix  B 

Doppler  Geometry 


The  calculation  of  the  plasma  velocity  component  radial  to  the  radar  involves 
the  solution  of  two  spherical  triangles  (see  Figure  Bl).  In  these  triangles 
M  =  geomagnetic  pole 
G  =  geographic  pole 

P  =  point  from  which  radar  energy  is  being  backseattered 
S  =  radar  site  location 
A  =  magnetic  declination 
B  =  magnetic  longitude  of  P 
<t> p  =  colatitude  of  P 
=  colatitude  of  S 
R  =  clutter  range 
a  =  clutter  azimuth 
a  =  magnetic  colatitude  of  G 
B  =  magnetic  colatitude  of  P. 

The  velocity  components  (V^  and  V^,)  are  specified  by  the  model  of  Sections  5.  1 
and  5.  2,  relative  to  M. 

Solving  triangle  MPG  yields  A, 

sin  B  sin  A  sin  0p 

tan  A  -  - c - - — 

cos  a  -  cos  b  cos  0p 

where  a,  b,  c  are  measured  in  angular  units  on  the  earth's  surface. 
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Figure  Bl.  Geometry  for  Doppler  Shift 
Calculations  in  Spherical  Coordinate 
System.  P  =  point  of  clutter  origination; 
S  =  radar  site;  G  =  geographic  pole; 

M  =  geomagnetic  pole 


Knowledge  of  A  permits  the  velocity  vector  to  be  transformed  into  the  geo¬ 
detic  system. 

Solving  triangle  SGP  yields  5,  the  angle  required  to  compute  the  velocity 
component  radial  to  the  radar  from  the  geodetic  components, 

sin  a  sin  9g  sin  R 
tan  6  '  (cos  4>s  -  cos  9 p  cos  R) 

where  R  is  measured  in  angular  units. 

The  velocity  vector  may  be  represented  in  the  geomagnetic  system  by 


V  =  (VVN+  VL2]A 


where 


Then  the  radial  velocity  component  relative  to  the  radar  (S)  is 
VR  =  V  cos  P 
where 

p  =  e  +  a  +  5  . 


Appendix  C 

Back  scatter  Volume 

Computation  of  backscatter  volume  in  the  presence  of  ionospheric  refraction 
(see  Figure  Cl). 

The  figure  depicts  two  rays  leaving  the  radar  (R)  at  an  angular  separation  (At ) 
in  the  vertical  plane,  and  undergoing  ionospheric  refraction.  The  magnetic  field 
vector  H  is  illustrated  at  three  adjacent  points  on  each  ray;  the  central  of  each  set 
of  three  points,  A  and  B  respectively  are  the  locations  of  geomagnetic  orthogonality 
top  =  90°)  corresponding  to  the  leading  and  trailing  edges  of  the  transmitted  pulse. 
Thus  the  mean  of  the  group  path  lengths  (along  the  ray  trajectories)  RA  and  RB 
correspond  to  the  radar  range  and  AB  defines  one  dimension  of  the  scattered  volume, 
determined  by  the  radar  resolution. 

The  scatter  volume  dimension  normal  to  AB,  in  the  vertical  plane  is  deter¬ 
mined  by  magnetic  aspect  sensitivity.  Let  C,  D,  and  E,  F  respectively  be  pairs 
of  points  on  either  side  of  the  orthogonality  points  at  which  the  backscatter  cross 
section  is  reduced  by  a  specific  amount  (say  6  dB)  from  its  maximum  value  at 
orthogonality.  CD  and  EF  are  then  given  by  2  •  Aip  where  P  is  group  path 

length  and  A\p  is  the  departure  of  the  geomagnetic  aspect  angle  from  orthogonality 
which  is  necessary  to  reduce  the  backscatter  cross  section  by  the  required  amount. 
A\p~±  l*.  for  6  dB  reduction. 

The  third  component  of  scatter  volume  is  the  azimuthal  component,  RA  0  where 

R  =  radar  range, 

Ap  =  azimuthal  antenna  beamwidth. 
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Figure  Cl.  Definition  of 
Scatter  Volume  for 
Refracted  Ionospheric 
Propagation 
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The  factor  F  in  this  expression  is  included  to  take  account  of  ionospheric  focussing. 

The  parameters  At,  and  F  must  be  determined  by  numerical  (or  in  some 
special  cases,  analytical)  ray  tracing  in  modelled  ionospheres. 

Figure  C2  shows  the  upper  part  of  a  computed  ray  trajectory  in  the  plane  of 
propagation  and  also  the  geomagnetic  aspect  angle  variation  along  the  trajectory. 

By  appropriate  scaling  from  this  figure,  the  value  of  ^  at  if/  =  90°  can  be  esti¬ 
mated  as  4.  0  km/degree.  Figure  C3  shows  the  variation  with  group  path  and  alti¬ 
tude  of  the  geomagnetic  orthogonality  point  (solid  line)  and  the  value  of  ^  (dashed 
line)  at  different  elevation  angles  of  take-off  at  the  transmitter.  This  gives  an  in¬ 
dication  of  the  amount  of  variability  to  be  expected  for  a  typical  case  (in  this  instance, 
6  MHz  propagation  from  a  site  at  61°  CGL  at  0°  azimuth  at  00  UT  in  January  at 
sunspot  number  100). 

Figure  C4  shows  the  group  path  to  the  geomagnetic  orthogonality  point  (for  the 
same  set  of  conditions  as  the  previous  two  figures)  as  a  function  of  elevation  angle 
at  the  radar.  Shown  also  is  the  range  of  elevation  angles  (At  =0.  9°)  corresponding 
to  a  250  fisec  pulse  at  a  radar  range  of  1000  km. 

The  focussing  factor,  F.  is  highly  variable  and  difficult  to  compute.  For  some 
representative  cases,  values  of  F ~  10  were  found.  The  significance  of  the  focussing 
factor  may  be  seen  by  noting  that  the  vertical  dimension  of  the  scattering  volume 
estimated  from  Figure  C3  is  much  less  than  RAt,  which  is  the  equivalent  vertical 
dimension  for  straight  line  radio  propagation. 
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Figure  C4.  Variation  Along  Hay  Trajectory  of  Group  Path  With  Elevation  Angle 
of  Take-off.  Showing  the  Kar.ge  of  Elevation  Angles,  at  1000  km  Kadar  Kange, 
Defined  by  the  Radar  Pulse 


Appendix  D 
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F-Region  Volumetric  Cron  Section 


The  probability  of  occurrence  of  F-region  irregularity,  p(F)  was  modelled  in 
Section  4.  2  in  terms  of  the  probability  of  spread-F  as  observed  by  vertical  inci¬ 
dence  ionosondes.  This  model  can  be  used  to  estimate  F-region  volume  reflectivity 
by  "calibrating"  it  with  Polar  Fox  II  cross  section  data.  Figures  16-19  show 
probability  distributions  of  F-region  "pseudo  volumetric"  cross  section  measured 
by  Polar  Fox  II  during  nighttime  hours  of  October  1972,  at  different  frequencies. 
Using  the  model  of  Section  4.  2,  the  probability  of  F-region  irregularity  in  the 
Polar  Fox  coverage  area  in  October  is  approximately  90  percent.  From  Figure  19 
the  value  of  10  MHz  cross  section  which  is  exceeded  90  percent  of  the  time  is 
-94  dB  m  / m  .  Thus  in  order  to  maintain  the  measured  difference  between  E-  and 

F-region  median  values  of  cross  section,  for  October  nighttime,  the  50  percent 

2  3 

value  of  p(F)  must  be  equated  to  -94  dB  m  /m  ,  for  average  Kp. 

The  procedure  adopted  here  makes  the  not  unreasonable  assumption  that  the 
volume  cross  section  is  proportional  to  p(F)  and  that  both  quantities  have  similar 
probability  distributions.  The  disadvantages  in  the  use  of  spread-F  incidence  as 
a  measure  of  backscatter  cross  section  is  that  it  is  too  crude  a  parameter,  having 
too  small  a  dynamic  range  and  too  easily  saturated.  Thus  the  range  of  estimated 
backscatter  cross  sections  may  be  somewhat  truncated. 
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Appendix  E 

Simple  Doppler  Shift  Models 

Section  5  of  this  report  contains  the  development  of  a  detailed  model  for  the 
Doppler  Shift  experienced  by  a  radar  signal  backscattered  from  auroral  irregulari¬ 
ties.  This  model  is  much  more  complex  than  previous  models  used  in  connection 
with  the  radar  aurora  and  may  in  some  applications  be  too  cumbersome  to  imple¬ 
ment  or  be  unnecessarily  precise.  Simpler  models  have  been  used  based  on  the 

assumption  that  radar  backscatter  originates  in  a  line  current  segment  (for  example. 
52 

Unger  et  al.  PARL  report).  These  models  are  based  on  the  common  assumption 
that  the  zonal  component  of  velocity  greatly  exceeds  the  meridional  component. 
Figure  El  shows  in  the  upper  diagram  a  straight  line  velocity  vector  (V),  showing 
the  reversal  of  the  radial  component  (V  )  as  the  azimuth  (a)  goes  through  zero  at 
true  north. 

V  =  V  sin  a  . 
r 

A  more  accurate  model,  but  one  still  having  the  advantage  of  simplicity  is 
illustrated  in  the  lower  diagram,  where  the  velocity  (V)  is  assumed  tangent  to  a 
line  of  constant  geomagnetic  latitude,  which  is  the  approximate  behavior  of  the 
auroral  electrojet. 

If 

4-  =  co-latitude  of  this  line, 
c 

4>g  =  co-latitude  of  radar  site,  0  , 
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Figure  E2  shows  the  azimuthal  variation  of  Doppler  Shift  derived  from  such  a 

model  with  4  =  22°,  4  =  29.  5°,  and  V  =  400  m/sec,  for  a  frequency  of  10  MHz. 
c  s 

Note  that  there  is  a  maximum  azimuth 


at  which  the  line-of- sight  becomes  tangent  to  the  velocity  vector  locus.  The 
tendency  of  the  Doppler  vs  azimuth  curve  to  approach  horizontal  in  this  model,  at 
large  azimuths,  is  frequently  observed  in  experimental  studies  of  auroral  back- 
scatter. 

The  straight-line  model  may  be  more  appropriate  for  clutter  deep  within  the 
polar  cap,  while  the  circular  model  is  more  appropriate  for  clutter  in  the  vicinity 
of  the  auroral  oval. 


Figure  E2.  Doppler  Shift  Calculated  on  Basis  of  Curved  Trajectory 
Model  in  Figure  El,  as  a  Function  of  Magnetic  Azimuth 
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Appendix  F 

Auroral  E-Layer  Effect 


Up  until  this  point  it  has  been  assumed  that  the  nighttime  ionospheric  refraction, 
below  the  altitude  of  E-region  irregularity  maximum,  is  small.  This  assumption 
permits  the  use  of  straight-line  geometrical  propagation  to  determine  the  geomag¬ 
netic  aspect  angle  at  100  km  altitude.  It  is  known,  however,  that  a  thick  refracting 
layer  of  ionization  is  often  present  on  the  night  side  of  the  auroral  oval.  The  peak 
electron  density  of  this  layer  varies  sporadically,  and  the  location  of  the  layer 

varies  with  geomagnetic  activity.  Recent  studies,  using  the  Chatanika  incoherent 
88 

scatter  radar  in  Alaska  have  revealed  details  of  this  auroral  E-layer  which  per¬ 
mit  construction  of  a  representative  model,  through  which  ray  tracing  can  be 
performed  to  establish  the  probable  regions  of  maximum  E-region  auroral  clutter. 

The  model  consists  of  a  torus-shaped  region  of  ionization  having  a  parabolic 
vertical  distribution  of  ionization  centered  at  110  km  altitude  with  10  km  semi¬ 
thickness.  The  latitudinal  effect  is  described  by  another  parabola. 

The  peak  density  of  the  layer  was  chosen  to  correspond  to  a  critical  frequency 
of  3  MHz,  which  is  supported  by  Chatanika  radar  measurements  for  low-to-moder- 
ate  magnetic  activity  levels.  Numerical  ray  tracing  was  performed  at  6  MHz  and 
12  Mhz  in  the  model  ionosphere  described  here,  using  the  Jensen-Cain  6th  order 
expansion  of  the  geomagnetic  field  to  determine  geomagnetic  aspect  angles.  The 

88.  Vondrak,  R.  R. ,  Smith,  G. ,  Hatfield,  V.  E. ,  Tsunoda,  R.  T. ,  Frank,  V.  R. , 

and  Perrault.  P.  D.  (1978)  Chatanika  Model  of  the  High- Latitude  Ionosphere, 
for  Application  to  HF  Propagation  Prediction,  KADC-Tk-78-7,  AbA053154. 
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models  were  constructed  for  =  0,  2,  4.  Figure  FI  shows  contours  of  90°  aspect 
angle  resulting  from  this  simulation.  The  dashed  lines  are  for  12  MHz  and  the 


solid  lines  for  6  Mhz  propagation.  These  contours  may  be  compared  to  the  90° 
contour  in  Figure  21,  which  is  for  unrefracted  propagation.  The  combination  of 
the  two  sets  of  contours  determines  the  total  effect  of  E-region  clutter. 


-  6  MHz 
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Figure  FI.  Contours  of  Calculated  Geomagnetic  Orthogonality  at  Altitude  110  km 
Using  an  Ionospheric  Model  Containing  an  Aurural  E-layer.  Dashed  lines  are  for 
12  MHz  and  solid  lines  for  6  MHz.  Site  is  located  at  61s  CGL 


MISSION 
of 

Rome  Air  Development  Center 

RAVC  plant  and  executes  research,  development,  test  and 
selected  acquisition  programs  In  support  of  Command,  Control 
Communications  and  Intelligence  (C3I)  activities.  Technical 
and  engineering  support  within  areas  of  technical  competence 
is  provided  to  ESV  Program  Offices  IPOs )  and  other  BSD 
elements.  The  principal  technical  mission  areas  are 
cormunicatf.ons,  electromagnetic  guidance  and  control,  sur¬ 
veillance  0$  ground  and  aerospace  objects,  Intelligence  data 
collection  and  handling ',  information  system  technology, 
ionospheric  propagation,  solid  state  sciences ,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


Printed  by 

United  States  Air  Force 
Hanscom  AFB,  Mass.  01731 


